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When holes are doped into an antiferromagnetic insulator they form a slowly fluctuating array of 
"topological defects" (metallic stripes) in which the motion of the holes exhibits a self-organized 
quasi one-dimensional electronic character. The accompanying lateral confinement of the interven- 
ing Mott-insulating regions induces a spin gap or pseudogap in the environment of the stripes. We 
present a theory of underdoped high temperature superconductors and show that there is a local 
separation of spin and charge, and that the mobile holes on an individual stripe acquire a spin gap 
via pair hopping between the stripe and its environment; i.e. via a magnetic analog of the usual 
superconducting proximity effect. In this way a high pairing scale without a large mass renormal- 
ization is established despite the strong Coulomb repulsion between the holes. Thus the mechanism 
of pairing is the generation of a spin gap in spatially-confined Mott-msulating regions of the mate- 
rial in the proximity of the metallic stripes. At non-vanishing stripe densities, Josephson coupling 
between stripes produces a dimensional crossover to a state with long-range superconducting phase 
coherence. This picture is established by obtaining exact and well-controlled approximate solutions 
of a model of a one-dimensional electron gas in an active environment. An extended discussion of 
the experimental evidence supporting the relevance of these results to the cuprate superconductors 
is given. 



I. INTRODUCTION 



Superconductivity in metals is the result of two dis- 
tinct quantum phenomena, pairing and long-range phase 
coherence. In conventional homogeneous superconduc- 
tors the phase stiffness is so great that these two phe- 
nomena occur simultaneously. On the other hand, in 
granular superconductors and Josephson junction arrays, 
pairing occurs at the bulk transition temperature of the 
constituent metal, while long-range phase coherence oc- 
curs, if at all, at a much lower temperature characteris- 
tic of the Josephson coupling between superconducting 
grains. High temperature superconductivity is hard 
to achieve, even in theory, because it requires that both 
scales be elevated simultaneously-yet they are usually in- 
compatible. Consider, for example, the strong-coupling 
limit of the negative U Hubbard model ||] or the Hol- 
stein model Pairs have a large binding energy but, 
typically, they Bose condense at a very low temperature 
because of the large effective mass of a tightly bound 
pair. (The effective mass is proportional to \U\ in the 
Hubbard model and is exponentially large in the Holstein 
model.) A similar issue arises if the strong pairing occurs 
at specific locations in the lattice (negative-f7 centers) ; in 
certain limits this problem may be mapped into a Kondo 
lattice Ql, which displays heavy-fermion behavior. 

A second problem for achieving high temperature su- 
perconductivity is that strong effective attractions, which 
might be expected to produce a high pairing scale, typ- 
ically lead to lattice instabilities, charge or spin density 
wave order, or two-phase (gas-liquid or phase separated) 
states B. Here the problem is that the system either 



becomes an insulator or, if it remains metallic, the resid- 
ual attraction is typically weak. In the neighborhood of 
such an ordered state there is a low-lying collective mode 
whose exchange is favorable for superconductivity, but 
the superconducting transition temperature is depressed 
by vertex corrections and also because the density of 
states may be reduced by the development of a pseudo- 
gap. 

A third (widely ignored) problem is how to achieve a 
high pairing scale at all in the presence of the repulsive 
Coulomb interaction, especially in a doped Mott insula- 
tor in which there is poor screening. A small coherence 
length (or pair size) implies that neither retardation, nor 
a long-range attractive interaction is effective in over- 
coming the bare Coulomb repulsion. Indeed, in the high 
temperature superconductors, angle resolved photoemis- 
sion spectroscopy (ARPES) suggests that the energy 
gap (and hence the pairing force) is a maximum for holes 
separated by one lattice spacing, where the bare Coulomb 
interaction is very large. 

In short, superconductivity typically occurs at low 
temperatures: if any attractive interaction is weak the 
pairing energy is small; if it is strong the coherence 
scale is suppressed or the system is otherwise unsta- 
ble. When this is coupled with the problem presented 
by the Coulomb force in a doped Mott insulator, the oc- 
currence of high temperature superconductivity in the 
cuprate perovskites is even more remarkable. Indeed, 
there is evidence [|-10| that these materials live in a re- 
gion of delicate balance between pairing and phase coher- 
ence: in "underdoped" and "optimally doped" materials, 
the onset of superconductivity is controlled by phase co- 
herence, and occurs well below the pairing temperature, 
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while in "overdoped" materials pairing and phase coher- 
ence take place at more or less the same temperature, as 
in more conventional superconductors. (See Fig. ^) If 
we accept this point of view, then we can approach the 
problem of understanding the mechanism of high tem- 
perature superconductivity from the underdoped side by 
addressing three separate questions: i) What gives rise to 
the large temperature scale for pairing, or in other words 
for superconductivity on a local scale? ii) How can the 
system avoid the detrimental effects of strong pairing on 
global phase coherence? {i.e. large mass renormaliza- 
tions.) iii) How can high temperature superconductivity 
with a short coherence length coexist with poor screening 
of the Coulomb interaction? 




FIG. 1. Theoretical sketch of the phase diagram for a 
high temperature superconductor in the doping-temperature 
plane. The solid lines represent phase transitions and the 
shaded areas crossovers. Tjv marks the transition to an an- 
tiferromagnetically ordered insulating state, and Tc the tran- 
sition to the superconducting state. Tj* marks the crossover 
temperature at which charge inhomogeneities (stripes) form 
and correspondingly local antiferromagnetic correlations de- 
velop in the insulating regions; the present paper is primarily 
concerned with the region between and somewhat above 
Tc, where the developing correlations are primarily confined 
to the neighborhood of an individual stripe. T2* marks the 
temperature scale at which a spin gap develops in the IDEG, 
and correspondingly the local superconducting susceptibihty 
begins to diverge. Here, T^, which is approximately 1/2 the 
antiferromagnetic exchange energy, marks the temperature at 
which the antiferromagnetic correlation length in the undoped 
antiferromagnet is equal to two or three lattice constants. For 
further discussion, especially concerning the experimental jus- 
tification for this figure, see Sec. IXC. 

Here we shall argue that the high temperature super- 
conductors resolve these problems in a unique manner: 1) 
The tendency of an antiferromagnet to expel holes ||l^] 



leads to the formation of hole-rich and hole-free regions 
Ip^ . For neutral holes this leads to a uniform instability 
(phase separation) jl^ but, for charged holes, the com- 
petition with the long-range part of the Coulomb inter- 
action generates a dynamical local charge inhomogcne- 
ity, in which the mobile holes are typically confined in 
"charged stripes", separated by elongated regions of in- 
sulating antiferromagnet [l^-|lj] . This self-organized col- 
lective structure, which we have named topological dop- 
ing is a general feature of doped Mott insulators, 
and it produces a locally quasi one-dimensional electronic 
character since, the electronic coupling between stripes 
falls exponentially with the distance between them [ p7[ . 
2) In a locally-striped structure, there is separation of 
spin and charge, as in the one-dimensional electron gas 
Ip^ (IDEG). Hence "pairing" is the formation of a spin 
gap, while the superfluid phase stiffness (i.e. the super- 
fluid density divided by the effective mass) is a property 
of the collective charge modes |19-2|]. 3) A large spin 
gap (or spin pseudogap) arises naturally in a spatially- 
confined, hole-free region, such as the medium between 
stripes. This effect is well documented for spin ladders 
p2t , and for spin chains with sufficient frustration |^,|^ . 
The important point is that the spin gap does not confiict 
with the Coulomb interaction since the energetic cost of 
having localized holes in Cu Zd orbitals has been paid 
in the formation of the material. 4) The spin degrees of 
freedom of the IDEG acquire a spin gap by pair hopping 
between the stripe and the antiferromagnetic environ- 
ment. (Single particle tunnelling is irrelevant |^^.) At 
the same time, because of the local separation of spin 
and charge, the spin-gap fixed point is stable even in the 
presence of strong Coulomb interactions, and there is no 
mass renormalization to depress the onset of phase co- 
herence, so the superconducting susceptibility diverges 
strongly below this temperature . 

In summary, the "mechanism" of high temperature su- 
perconductivity is a form of magnetic proximity effect 
in which a spin gap is generated in Mott-insulating an- 
tiferromagnetic regions through spatial confinement by 
charge stripes, and communicated to the stripes by pair 
hopping. The mobile holes on the stripes have the large 
phase stiffness required for a high superconducting tran- 
sition temperature. 

The formation of a spin gap in the IDEG may be re- 
garded as a pairing of "spinous", i.e. the neutral, spin- 
1/2 soliton excitations which occur in the low energy 
spectrum of the IDEG and a number of one-dimensional 
quantum antiferromagnets. Indeed, local inhomogeneity 
provides a realization of some of the earlier ideas p^ ] 
involving spin-charge separation in the high tempera- 
ture superconductors and the concept of a spin liquid, 
by which we mean a quantum disordered system (i.e. 
with unbroken spin-rotation symmetry) which supports 
spinous in its physical spectrum. However, we emphasise 
that previous ideas relied on a putative two-dimensional 
spin-liquid fixed point, while here we are dealing with a 
locally one-dimensional system, for which it is well es- 
tablished 1^,0 that separation of spin and charge |l^ ] 
occurs generically, and there exists a "paired spin-liquid" 
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phase, i.e. a spin-liquid with a finite gap or pseudogap in 
the spinon spectrum. (See discussion in Appendix C.) In 
the strictest sense then, both are intermediate-distance 
effects [0 which occur below a dimensional-crossover 
scale to two (or three) dimensional physics. 

We thus view the emergence of high temperature su- 
perconductivity as a three-stage process, which can be 
described in renormalization group language in terms of 
the influence of three fixed points. At high temperatures, 
the "avoided critical phenomena" |^ associated with 
frustrated phase separation, govern the emergence of 
the self-organized, quasi one-dimensional structures. At 
intermediate temperatures, the one-dimensional paired 
spin liquid fixed point controls the pairing scale, and the 
growth of local superconducting (and CDW) correlations. 
Finally, at low temperatures, a two (or three) dimen- 
sional fixed point determines the long-distance physics 
and the ultimate superconducting or insulating behavior 
of the system. 

Our proposed mechanism implies the existence of two 
crossover scales above Tc in underdoped materials, as 
shown in Fig. 1: a high temperature scale, at which 
local stripe order and antiferromagnetic correlations de- 
velop, and a lower temperature at which local pairing 
(spin gap) and significant superconducting correlations 
appear on individual charge stripes. Tc itself, is then de- 
termined by the Josephson coupling between stripes, i.e. 
by the onset of global phase coherence Q . 

The local charge inhomogeneity which is a central 
feature of our model has substantial support from ex- 
periment. In the past few years charge ordering has 
been discovered in a number of layered oxides, such 
as La2-KSra;Ni04+<5 Q and Lao.5Sri.5Mn04 and 
there is considerable experimental evidence showing that 
the high temperature superconductors display a coex- 
istence of superconductivity and charge inhomogeneity. 
In particular, the efficient destruction of the antiferro- 
magnetic order | |40[ | of the parent insulating state is a 
consequence of topological doping jl^ , in which the mo- 
bile holes form metallic stripes that are antiphase domain 
walls for the spins. The stripes may be ordered ||4l[| (as in 
Lai.6-2;Ndo.4Sr^Cu04), dynamically fluctuating [ p]j42| 
(as in optimally-doped La2-a;Sra;Cu04), or pinned and 
meandering |43| (as in lightly doped La2-2:Sra;Cu04). 
Thus, we consider the existence of local metallic stripes 
(at least in the La2Cu04 family of high temperature su- 
perconductors) to be an experimental fact. Evidence of 
specific charge fiuctuations in any family of cuprate su- 
perconductors suggests that they are an important in- 
gredient in the theory of high temperature superconduc- 
tivity. However neutron scattering data also sug- 
gest that there are similar, but more disordered, struc- 
tures in underdoped YBa2Cu307_5. An analysis 
of ARPES experiments on Bi2Sr2CaCu208+5 leads to a 
similar conclusion p5[ . 

The systematics of phase fluctuations Q, mentioned 
above, strongly suggests that pairing on a high energy 
scale does not require interaction between metallic charge 
stripes, although Tc is certainly controlled by the Joseph- 
son coupling required to establish phase coherence for an 



array of stripes. Consequently, it should be possible to 
understand the mechanism of pairing from the behavior 
of a single stripe, modelled as a IDEG coupled to the 
various low- lying states of an insulating environment. A 
complete discussion of this problem is a substantial gen- 
eralization of the theory of the one-dimensional electron 
gas p6[ which will be considered more completely in a 
subsequent publication Here it will be shown that, 
for the high temperature superconductors, the most im- 
portant process is the hopping of a pair of holes from 
the stripe into the antiferromagnetic environment, which 
also may be regarded as a coherent form of transverse 
stripe fluctuation. It will be shown that the stripe de- 
velops a spin gap which, in this model, corresponds to 
pairing without phase coherence. We consider two situ- 
ations: a) the antiferromagnetic environment has a pre- 
existing spin gap or spin pseudogap because of its finite 
spatial dimensions and b) pair hopping produces a 
spin gap in both the stripe and the environment. In the 
first case, we find that an induced spin gap in the IDEG 
and the consequent divergent superconducting fluctua- 
tions are a robust consequence of the coupling to the en- 
vironment. The second case requires a sufficiently strong 
(and possibly unphysical) Coulomb interaction between 
holes on the stripe and holes in the environment for pair 
tunnelling to be relevant. 

Although the existence of two distinct regions, the 
stripe and the antiferromagnetic environment, provides 
a potential escape from some of the limitations on the 
superconducting transition temperature Tc, it is not a 
priori obvious that a large mass renormalization can be 
avoided. Indeed, the model we shall study is closely 
related to Kondo lattice models 0], for which heavy- 
fermion behavior or large mass renormalization is the 
primary consequence of the strong interactions. How- 
ever we find that, for stripes in an antiferromagnet (as for 
one- dimensional Kondo and orbital Kondo lattice models 
1^,^), the analog of heavy- fermion physics is reflected 
solely in the the spin degrees of freedom while for the 
charge modes, and hence the superfluid phase stiffness, 
the mass is not renormalized! 

In some respects, what we are doing is analogous to 
working out the renormalization of the electron self en- 
ergy by the coupling to phonons. However, the calcula- 
tion is more complicated because, here, the elementary 
objects are strings of charge (stripes) in a polarizable 
medium that profoundly influences their internal struc- 
ture. Fluctuating stripes are of flnite length but the solu- 
tion of the inflnite IDEG may be used if they are longer 
than the spin gap length scale, which is a few lattice 
spacings Q. 

Of course, at higher hole concentrations, the calcu- 
lation must be modified to take account of the inter- 
action between the stripes, especially to obtain long- 
range superconducting order. In general terms, it is 
fairly straightforward to see how global superconductiv- 
ity arises in a system with a small but finite density of 
ordered or slowly- fluctuating stripes, as found in under- 
doped members of the La2-xSra;Cu04 family of super- 
conductors. Indeed, an analysis of neutron scattering 
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and thermodynamic data for underdoped and optimally 
doped La2-£!;Sra;Cu04 j4|] suggests that Tc is propor- 
tional to the product of the Drude weight of the holes on 
a stripe and the stripe concentration c^. 

An interesting feature of our model is the interplay be- 
tween the short-distance physics associated with the fluc- 
tuating stripes and the ultimate long-range order that is 
established in a given material. We shall show that both 
superconducting and charge density wave correlations de- 
velop on a given stripe. However, they compete at longer 
length scales, although they may coexist in certain re- 
gions of the phase diagram. Also it follows from general 
principles that, locally, the singlet superconducting order 
parameter will be a strong admixture of extended-s and 
Ax2_y2 states. Ultimately, in tetragonal materials, the 
order parameter must have a pure symmetry, but the 
way in which it emerges from the short-distance physics 
is very different from more conventional routes. 

This paper is quite long and, in parts, rather techni- 
cal. It addresses the purely theoretical problem of con- 
structing and solving a general model of a IDEG in an 
active environment. At the same time, we wish to re- 
port progress on the key problem of understanding the 
mechansim of high temperature superconductivity in the 
cuprate superconductors. To compensate, we have at- 
tempted to make the various sections as self-contained as 
possible, and to indicate which sections can be skipped by 
the reader with a more focussed interest in the problem. 

A rather general model of the interactiong IDEG in an 
active environment is introduced in Sec. II. The model 
is bosonized in Sec. Ill, and various formal transforma- 
tions that are useful for later analysis are described; this 
section also contains a discussion of which of the allowed 
interactions in the model are unimportant for our pur- 
poses, and so can be ignored. In Sec. IV, we define a 
simplified "pseudospin" model of the charge excitations 
of the environment, and argue that it exhibits the same 
low-energy physics as the general model. Sec. V con- 
tains a discussion of exact results for the zero temper- 
ature properties of the pseudospin model, which among 
other things exhibits the spin-gap proximity effect, and 
the generation of a paired spin liquid state of the IDEG, 
even in the presence of arbitrarily-strong forward scat- 
tering. Section VI reports the results of a controlled ap- 
proximate solution of the pseudospin model for a wide 
range of temperatures and coupling constants; in partic- 
ular, various crossover temperatures to spin-gap behavior 
are identified, and their dependence on the interactions 
in the model are determined. In Sec. VII, we return 
to the problem of the charge degrees of freedom of the 
IDEG, and consider the effects of umklapp scattering in 
conditions of near commensurability, and the effects of 
an externally applied potential. In Sec. VIII, we digress 
slightly to consider the effects of a "spin-gap center" on 
the local properties of a Fermi liquid. Finally, in Sec. 
IX, we summarize our results and discuss experimental 
implications and predictions for the high temperature su- 
perconductors. In this section, we also suggest some nu- 
merical calculations to test the major ideas. The reader 
who is primarily interested in a discussion of results may 



skip directly to Sec. IX. In addition. Appendix A recasts 
some of the present discussion in the familiar language 
of the perturbative renormalization group for the IDEG, 
Appendix B contains an analysis of the symmetries of 
the model, and an explicit construction of the non-local 
order parameter which characterizes "local pairing" , and 
Appendix C discusses the precise nature of the paired- 
spin-liquid state, and gives concrete examples of model 
systems which exhibit this state. 



II. THE IDEG IN AN ACTIVE ENVIRONMENT 
A. The problem and the solution strategy 

It has long been realized that the low energy properties 
of a one dimensional electron gas (IDEG), and indeed of 
a wide variety of other interacting one dimensional sys- 
tems, are equivalent to those of the simplest field theory 
of interacting electrons, characterized by a small number 
of potentially relevant interactions between electrons at 
the Fermi surface. In this section we address the problem 
of a IDEG in an "active" environment, one that possesses 
its own low-energy excitations which couple to the IDEG, 
but is insulating so that the electrons of the IDEG may 
make excursions into the environment, but ultimately 
return. The environment in which we are interested is 
antiferromagnetic, so it may have low-energy spin exci- 
tations. It will also have low-energy charge excitations in 
which holes make excursions from the metallic stripe into 
the environment. Their energy is low because frustrated 
phase separation, which generates metallic stripes in the 
first place, involves a delicate balance of Coulomb and 
magnetic energies. 

This problem can be addressed in several distinct 
ways. In the present paper, we make extensive use of a 
renormalization group strategy involving exact solutions 
of solvable models, together with a sophisticated approx- 
imate calculation, in which the fluctuations of the IDEG 
and the environment are solved exactly, but the coupling 
between them is treated in a mean-field approximation. 
We also give physical estimates of the values of the vari- 
ous coupling constants that enter the model, and present 
strong physical arguments to show that the physical sys- 
tems of interest will lie in the "basin of attraction" of the 
strong-coupling fixed point that governs the behavior of 
the solvable models. In Section IX, we will also out- 
line some simple one-dimensional lattice models which 
are amenable to numerical solution, and are expected to 
exhibit the mechanism described in this paper. 



B. The general model 

To begin with, we consider a very general model of a 
IDEG coupled to an environment. The initial form of 
the model is microscopically realistic. It will be assumed 
that the environment itself is a one dimensional system 
with a charge gap (since it is an insulating matrix) which 
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may or may not have a spin gap. We thus consider the 
Hamihonian to be of the form 

/C30 
dx [HiDEG + 'Henv + 'Hint + Ti-coul] ■ (1) 
-oo 

The bare Hamiltonian density of the IDEG is 

Ti-iDEG = H-o + Hi. (2) 

Here Hq is the Hamiltonian of a non-interacting IDEG, 
which in the continuum hmit can be written (with h = 1) 

as 

a 

-IJ'^['^Laix)i^a.aix)] (3) 

where ipl^ ^ {x) creates an electron with z component of 
spin (T on the right or left moving branch of the Fermi 
surface for a = 1 or 2 respectively. Here, we have made 
a Galilean transformations to shift the Fermi points to 
fc = 0; factors involving the Fermi wave vector kp will be 
shown explicitly. Hi incorporates the electron-electron 
interactions within the IDEG and has the continuum 
form 

Hi^ V'lffV'L' "02,^7' V'i,<T 

+91 ■0L<tV'2,^'V'1,o-' V'2,a 
a,cr' 

+93 K,^V'la^2,iV^2,Te^'''""''^" + H.c] . (4) 

Here G is a reciprocal lattice vector and 93 is the cou- 
pling constant for umklapp scattering. When the IDEG 
is incommensurate {4:kp =/= G), the rapid phase oscilla- 
tions in the term proportional to 93 render it irrelevant 
in the renormalization group sense. However, near to 
commensurability, this term is responsible for the fact 
that the Drude weight is proportional to the density of 
doped holes, as we shall see. Typically, it will be as- 
sumed that the interactions are repulsive {91,92,93 > 0) 
although they may undergo significant renormalization 
by the coupling of the IDEG to the high energy excita- 
tions of the antiferromagnetic environment (which we do 
not consider explicitly). The parameters that describe 
the IDEG are thus the Fermi velocity, vf, the chem- 
ical potential, /x, the three coupling constants gi, and 
the "incommensurability", Akp — G. It should be em- 
phasised that this is a very general representation of the 
low-energy physics of a stripe in a CUO2 plane, and all 
details of the original microscopic model are contained in 
the values of the coupling constants gi. 

We have in mind the low-density limit of a stripe phase 
in which the Coulomb interaction on a given stripe is 
screened by the motion of charge on neighboring stripes, 
and so does not make a singular contribution to the for- 
ward scattering interaction, g2. Thus, for the time being, 
we will neglect the term Hcoui , although it will ultimately 



play a role in the dynamics of the superconducting phase 
I- 

Because the physics of interacting systems in one di- 
mension is ultimately so constrained, it is possible to 
model the Hamiltonian density of the environment as 
a second (distinct) interacting one dimensional electron 
gas. The Hamiltonian Henv has the same form as in Eqs. 
(^) and (^, except that fields and parameters will be 
marked with a super-tilde. However there are several im- 
portant differences in the parameters of the Hamiltonian: 
1) The environment is a Mott insulator. Consequently 
there is a strong commensurability energy (4^:^ = G and 
g3 is large), which produces a gap in the the charge de- 
grees of freedom of the environment. This also implies 
that kp is different from kp. 2) Because of the frustra- 
tion of the motion of holes in an antiferromagnet, po| the 
propagation velocity Vc for charge excitations in the envi- 
ronment is much smaller than the corresponding velocity 
in the IDEG. This is the primary manner in which the 
driving force for phase separation and stripe forma- 
tion appears in the model. 3) We shall consider 
three possibilities for the spin degrees of freedom of the 
environment, one in which there are gapless magnon-like 
excitations, and two in which there is a spin gap: a) The 
gapless state is realized by considering the model with 

> 0, in which case the environmental spin excitations 
are those of an antiferromagnetic spin-1/2 Heisenberg 
chain, b) A spin gap can occur with an accompanying 
spontaneous breaking of translational (chiral) symmetry 
(See Appendix B), which is realized by simply taking 

< 0, in which case the environmental spin excitations 
are those of a spin-1/2 Heisenberg chain with compet- 
ing nearest and next-nearest neighbor antiferromagnetic 
interactions, e.g. the Majumdar- Ghosh model. j2^] c) 
A spin-gap can occur without any accompanying broken 
symmetry, in the manner of the antiferromagnetic two 
leg, spin 1/2 Heisenberg ladder |2^; to model this sys- 
tem, we need to add a backscattering term to the environ- 
mental Hamiltonian (of the same form as in Eq. (|8^), 
below), although a better description can be attained in 
the bosonized form of the Hamiltonian, as discussed be- 
low. For our purposes, there is no significant difference in 
the implications of the two types of environmental spin 
gap, so for simplicity, we will perform our calculations for 
the case in which the spin gap is induced by a negative gi , 
and will use language to describe the physics that (prop- 
erly) does not distinguish the two types of environmental 
spin gap. 

Using well known results for the IDEG, it is possible 
to express these coupling constants in terms of the phys- 
ical variables which define the excitation spectrum of the 
environment: the spin and charge velocities, Vs and Vc, 
the charge gap Ac and the spin gap (if one exists) As, 
and the charge and spin correlation exponents (defined 
below), Kc and Kg. Since the environment is an insula- 
tor, we will always assume that Ac is large. We also must 
include the energy e to transfer charge from the IDEG 
to the environment. For the case of "p-type" doping, in 
which fi lies in the lower half of the environmental gap. 
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e/2 = Ac — [/I — /z] is the bare energy required to remove 
a quantum of charge from the environment and add it to 
the IDEG. We will be interested in the case < e <C Ac. 

Finally, we consider the coupling between the IDEG 
and the environment, for which spin-rotational invari- 
ance and conservation of momentum along the stripe di- 
rection severely limit the number of possible relevant in- 
teractions. Since the Fermi wave vector of the IDEG is 
incommensurate with the wave vector of any low energy 
excitation of the environment, we can neglect, as irrel- 
evant, terms which transfer momentum ±kp or ±2kp 
between the IDEG and the environment. For example 
there are no low energy single-particle hopping processes, 
even though, at the microscopic level, one might expect 
them to have the largest coupling term. Such processes 
are included implicitly as virtual intermediate states in 
constructing the effective low energy Hamiltonian. (We 
will return to this point briefly in the following section.) 
With this in mind, the most general form of the interac- 
tion Hamiltonian density, i.e. which keeps all potentially 
relevant terms, is 



Ti-int = J sis ■ js + VsS ■ S 

+ Jc.jc.jc + Vcpp 

air 5 



(5) 



where the small momentum transfer couplings involve 
the long-wavelength density fluctuations relative to the 
background charge density po 

p{x) - po = X! ['/'i,<tV'1,<t + V'l,aV'2,<T] , (6) 



the bare charge-current operator 

jc{.x) = ^ [V'UV'1,<T - ^l,a'^2,a\ , 



(7) 



the long-wavelength spin density operator 

(T.a' 

and the bare spin-current operator 

(7,(T' 

The corresponding operators for the environment are de- 
fined by the same equations, except that all quantities 
have a super-tilde. Note that we have chosen to express 
Ti-int in terms of the charge and spin current operators 
for the noninteracting system. The other contribution to 
Ti-int is the pair transfer terms 



Tipa,r=tsp[P^P + 'a.C.] 
1 



(10) 



where for the IDEG, is the usual singlet pair creation 
operator. 



PHx) = ^ [^l^{x)4jx)+4^^{x)^\jx)\ , (11) 

and Pm are the componenets of the triplet pair creation 
operator. 



Pi (x) = 



Po{x) = ^1^ {x)iPI^{x) - 2pl^{x)i^l^{x) 



P 



t - „/,t 



(12) 



III. BOSONIZATION OF THE MODEL 

In dealing with the problem of the IDEG in an ac- 
tive environment, it is useful to rewrite the model using 
the standard boson representation of Fermi fields in one 
dimension MM: 



1 



-.exp{i^x^a{x)} 



(13) 



where $A,a = 0F[6'^(x) ± (f>aix)] with "-" and 
corresponding to A = 1 and 2 respectively, 9c,{x) — 
J^^dx'Ila{x'), and <j)cr{x) and Ila{x) are canonically 
conjugate Bose fields, so that [0o-(a;), Ho-' (x')] = iS{x — 
x'). (0 and are thus dual to each other in the usual 
statistical mechanical sense of order and disorder vari- 
ables.) To take advantage of the separation of spin and 
charge the Hamiltonian will be expressed in terms 
of a spin field, ^^(a;) — [(j)] — (/)x]/V2, and a charge field, 
(j)c{x) = [(f)^^ + <f)i\/^/2 and their conjugate momenta: 
Tis{x) = Pi - Hj/y2 and Hc(a;) = [H| + HJ/V2. The 
charge and spin density and current operators may be 
written: 



p{x) = -\J^d^d)c (14) 

jcix) = \l\ He 

S^{x) 
S^{x) 

fsix) 

jf{x)^—exp{±iV2^es)sm[V2^(f>s]. (15) 



— A / — dx 
2tt 



exp(±iv27r0s) cos[v27r(/ 



27r 



In terms of these variables, the Hamiltonians of the 
stripe, the environment, and the small-momentum trans- 
fer coupling between the two may be written as a sum 
of a charge-only part and a spin-only part. However, the 
pair hopping terms Hpair introduces a coupling between 
spin and charge. Thus the total Hamiltonian may be 
written 



Ti — Tic + Tis + Tipair- 
We now consider the various contributions in turn. 



(16) 
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A. Spin Degrees of Preedom 



Ale = /tg + rtg + rtg, 



(21) 



The general form of the spin Hamiltonian is 

ni — 'T_/0 I a,/! I 



where 



(17) 



Here 



2J ~ 2V 

nl = — -IIs{x)ns{x) H -dj;(t)s{x)dx(j)s{x), 



and 



ni 



2.91 



(27ra)2 



(18) 



(19) 



(20) 



+._J|_cos[y8^^.] 



cos [V27r(6's — ^s)] cos [V27r(/)s] cos ['\/27r^s] 



+ ^ cos [V27r{es ~ 9s)] sin [V2^0.] sin [\/2^0,] 

Here Vg is the spin-wave velocity and Kg is the critical 
exponent |^l[] that specifies the location on a line of fixed 
points. Also Vs is given by Vg = 2vfKs/{K^ + 1). In the 
absence of coupling between the stripe and the environ- 
ment, the Hamiltonian is known to be correct for weak or 
strong coupling and for different forms of short-distance 
or high-energy cutoff [Q, although it may be necessary 
to perform some form of global renormalization to deter- 
mine Kg from the parameters of the initial Hamiltonian. 
For weak coupling, Kg is related to the bare Fermi veloc- 
ity vp and coupling constants as Kg = \f^^^^^- ^^o^' 

repulsive interactions (i.e. gi > 0) one finds Kg > 1. 

For the case in which gi is negative and relevant, in the 
renormalization group sense, there is a two fold degen- 
erate ground state, corresponding to the classical values 
= and 4>s = ^7r/2. (See Appendix B.) To rep- 
resent the case in which there is an environmental spin 
gap without symmetry breaking, we should add a term 
proportional to cos[-\/27r0s], which arises in a microscopic 
system with two spins per unit cell, such as a two-leg lad- 
der Is^ ]. This term (which may be generalized to allow 
any even number of spins per unit cell) is always relevant 
for repulsive interactions, so it always leads to a spin gap. 
As we shall see shortly, the important point is that a spin 
gap of whatever origin implies a quenching of the fluctu- 
ations of 04,. For a caveat on commensurability effects, 
see Sec. VH. 



B. Charge Degrees of Freedom 

The general form of the charge Hamiltonian is 



H^ = |[ifeH^ + -i-(a.0e)^] 

+ |[i^,n? + ^(a.^,)2] 



2Jr ~ 2Vr 
TT TT 



(22) 



(23) 



and 



= 



253 



(2™)2 

2g3 

'(27ra)2 



cos 



in - (4fcF - G)x] 



cos [-s/Stt^c] ^ ~ 9- 



(24) 



Here Vc is the charge velocity and Kc is the Luttinger 
liquid exponent |0, with Vc — 2vfKc/{K'^ + 1). For 
weak coupling, Kc is related to the bare Fermi velocity 

VF and coupling constants as Kc = y^ii^^^f^ i where 
gc — 91^ 2g2- For repulsive interactions < Kc < 1 {i.e. 
9c < 0). 



C. Spin-Charge Coupling 

Pair hopping between the stripe and the environment, 
as given by TLpair, destroys the separation of spin and 
charge and is the driving force for much of the interesting 
physics. Its bosonized form is given by 



H 



tg 



pair — 1 on 



Hp 



cos[\/27r( 



c)] cos[v 27r(/)s] cos[v27r(/>s] 



cos[\/2^(6ic - Be)] cos[\/2^(6is - Og)] 



■ cos[V27r(6'c - Be)] sin[\/27r0s] sin[\/27r0s] \ . (25) 



D. Which Terms Are Unimportant? 

The general model has numerous coupling constants, 
and so, for much of this paper, we focus on the terms that 
are most important for our purposes, and set the others 
to zero. Specifically, we drop those terms which are, in 
the renormalization group sense, irrelevant at the paired 
spin liquid fixed point. This argument simply shows that 
dropping these terms is self consistent. However, given 
the nature of the antiferromagnetic environment, there 
are strong arguments to show that these terms also are 
physically irrelevant, i.e. that the physical system lies 
in the basin of attraction of the paired spin liquid fixed 
point. 

To begin with, we examine the magnetic interactions, 
Jg and Vg in "Hint', these terms represent the interaction 
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between the ferromagnetic fluctuations in the two sub- 
systems. Since we are primarily interested in antiferro- 
magnetic systems, we do not expect these terms ever to 
be important. Of course, in the paired spin Uquid state, 
or more generahy in the presence of any sort of envi- 
ronmental spin gap, this can be seen directly from their 
dependence on 9s, which means that the corresponding 
correlation functions decay exponentially with distance 
or time, and are thus trivially irrelevant. The triplet 
pair-tunnelling term similarly depends on 9s, and cor- 
respondingly triplet pairing is generally expected to be 
important only in nearly ferromagnetic systems. There- 
fore, on both clear physical, and formal renormalization 
group grounds, it is safe to simplify our further discussion 
by taking 



Js^Vs 



tp 



0, 



(26) 



unless explicitly stated otherwise. 

Thus, in the case where there is strong incommen- 
surability between the values of kp in the two subsys- 
tems, and neither has significant ferromagnetic fluctua- 
tions, the only important interactions between the f DEG 
and the environment are tsp, Vc, and Jc- 

Away from half flUing, the renormalization of the umk- 
lapp scattering co upling constant 53 is cut off by the in- 
commensurability ||5^,Q, and for some purposes it may 
be dropped. However, this does not mean that umk- 
lapp scattering is unimportant for the low energy physics. 
Doping of holes into the Mott insulating state in one di- 
mension creates soliton excitations [^3|j5J] in the charge 
density with a mass governed by 53. There is a "doped- 
insulator" region in which these excitations control the 
Drude weight and the superfluid phase stiffness. In our 
stripe model of the cuprates, high temperature supercon- 
ductivity may occur within this region of doping. 

Finally, we address the non-linear term proportional 
to gi in Tij in Eq. (|20|). For repulsive interactions, i.e. 
for Ks > 1, this term is pcrturbatively irrelevant, and the 
renormalization group flows go to the flxed point gi = 
and Ks = 1. (See Appendix A.) Thus, so long as the 
bare interactions in the IDEG are not too large, it is 
reasonable to use the fixed point values 

51 = and Ks = 1 (27) 

for the effective low energy theory. 



E. Unitary Transformation 

We now introduce a unitary transformation which will 
be used in a number of ways to simplify the problem. 
The operator 



U\ = cxp[-iA / dx 9c{x)dx4>c{x)] 



(28) 



has the effect of shifting the fields 



Ulllc{x)Ux = U,{x) + XU,{x) 

Ulu,{x)Ux = n,(x) 
UlMx)Ux = Mx) - XMx)- 



(29) 



This transformation modifies the various charge interac- 
tions 



2 Kr 



Vc 

Jc AJc = Jc- ^XvcK, (30) 
and the velocities and exponent parameters 



Vc 

Kc 

Vc 

Kc 



Vcl 

Kcl 

Vcl 

Kc/l 



(31) 



where 



7 = 



X^icKc 
VcKc 



4AJe 
TTVcKc 



X^VcKc AXVcKc 



VcKc 



TTVc 



(32) 



1 . Perturbative Relevance of Pair Hopping 

The transformation ( p8| ) diagonalizes the quadratic 
part of the charge Hamiltonian + Hi provided 

. 2VcKr 



Jc = — 



TTVc 

VcKcKcVc 



Vc 



(33) 



We are now in a position to discuss the perturbative 
relevancy of pair hopping, which is the process that will 
generate a spin gap along the stripe. Here we have in 
mind the initial stage of renormalization, in which de- 
grees of freedom with energies large compared to the 
charge transfer energy, e, are eliminated. Thus it is rea- 
sonable to determine perturbative relevance relative to 
the quadratic piece of the Hamiltonian (See also 

Appendix A.) However other relevant pertubations, such 
as ^3, are important for the later stages of renormaliza- 
tion. Substitution of Eqs. ( ^3| ) into Eqs. ( |3^ ) gives 



41/2 v^KcKc 



4^2 KcKc 



1/2 



VrVc 



1/2 



(34) 
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Then the singletpair hopping operator Tipair is pertuba- 
tively relevant |46| if the exponent 



as 



satisfies a 



{l-Xf 



+ K, 



(35) 



sp 



J<c ^ 

< 1, and perturbatively irrelevant other- 



wise. Despite appearances, asp shares the property of 
the Hamiltonian that it is symmetric under interchange 
of Kc and K^. when = ^c- If interactions in the 
original model were set equal to zero, then all of the 
K's and 7's would be equal to 1, so that asp — 1, and 
pair hopping would be marginal. Repulsive interactions 
within the stripe and the environment increase the value 
of asp, since they make Ks,Ks > 1 and Kc,Kc < 1- 
This is physically reasonable because repulsive interac- 
tions within the stripe and the environment are unfavor- 
able for pairing. 

There are three effects which enhance the perturbative 
relevance of singlet pair hopping. First of all, it can be 
seen from Eqs. (|34| ) and (|35|), that a repulsive Vc de- 
creases the value of asp. Physically, this occurs because 
the charge density in the environment decreases in the 
vicinity of a pair in the IDEG; thus it is easier for the 
pair to hop. This effect is surely an important piece of the 
physics of pair hopping and it provides a way in which 
the Coulomb repulsion is favorable for pairing. But it 
cannot be the sole reason for the relevancy of singlet pair 
hopping unless Vc is greater than a suitable average of \gc\ 
and \gc\- As discussed in Appendix A, this can happen, 
in principle, if the character of the screening is just right, 
but it seems to be an insufficiently robust mechanism for 
a high temperature scale for pairing. 

Secondly, the frustration of the motion of holes in an 
antiferromagnet implies that the bare Fermi velocity vp 
of the environment is small, and hence Vc is small, which 
depresses the value of 7 (Eq. (|^) and the first contri- 
bution to asp in Eq. (35). 

Thirdly, if the environment has a preexisting spin gap, 
then one should set Kg — in the expression for asp; this 
substitution makes singlet pair hopping perturbatively 
relevant (i.e. asp < 1) for a wide range of the other 
parameters. A slightly weaker form of this route occurs 
if the environment has a spin pseudogap. For example 
it might have several gapped spin excitations and one 
gapless spin excitation, as in odd-leg ladders Then 
the i^s-term in asp should have a coefficient Wg < 1 equal 
to the weight of the gapless excitation in the pair hopping 
process. The elimination or reduction of Kg in I^q- ( |35| ) 
is the perturbative renormalization group manifestation 
of the proximity effect. 

It is important to note that transverse fluctuations of 
the stripe, together with the Coulomb interaction be- 
tween holes on the stripe and in the environment, in- 
crease the value of the superexchange coupling along 
neighboring bonds perpendicular to the stripe 
Clearly these processes decrease the value of Wg and are 
almost as effective as a full environmental spin gap for 
making pair hopping perturbatively relevant. Moreover 
the environment will vary along the length of a fluctuat- 
ing stripe, and singlet pair hopping may be relevant at 



some stripe locations ( "spin-gap centers" ) and irrelevant 
at others, where it may be neglected. This sort of local 
fluctuation is readily included in the pseudospin model 
introduced in the next section. 

The spin gap proximity effect, enhanced by a large Vc 
and small vp, gives a robust mechanism for the perturba- 
tive relevance of pair hopping for a wide and physically 
reasonable range of interactions. Similar conclusions can 
be drawn from examining the perturbative expression 
for the beta function for tgp in powers of the interaction 
strength, as is discussed in Appendix A. 



2. Composite Order Parameter 

In the rest of this paper, we shall use the canoni- 
cal transformation (|2^) in a slightly different way by 
taking A = 1, which is similar to the transformations 
employed [^,^ in the analysis of Kondo and orbital- 
Kondo arrays in one dimension. The special values of 
the coupling constants Vc and Jc for which the quadratic 
part of the charge Hamiltonian Tig is diagonalized at the 
point A = 1 are the analog of the Toulouse limit of the 
Kondo problem and the various decoupling lines of the 
multi-channel Kondo problem, and Kondo lattice prob- 
lems mill-Ill. 

For A = 1, the transformation eliminates the 9c depen- 



dence of UlTipairUi since Ul[9c — Oc]Ui = Oc- Remark- 
ably, this also implies that the transformed 9c is gauge 
invariant. Consequently it is possible to define a compos- 
ite superconducting order parameter |59 6^ in terms of 

Ui as, OcoTnp = t^iV'i,T'/'2,iC^i = {2TTay^ cxp[~iy/2TT{9c~ 
Oc + iips)], which can exhibit long-range order at zero 
temperature, despite the constraints of the Hohenberg- 
Coleman- Mermin- Wagner theorem for a conventional or- 
der parameter. Indeed, as discussed in Appendix B, long- 
range composite order implies a broken Z(2) symmetry, 
which, for lack of a better name, we call r symmetry. 

The transformation introduces a ipc dependence into 
the (73 term of 7i^, which complicates the analysis some- 
what although, as we shall see, it can be handled. How- 
ever, whenever 173 can be neglected, the unitary trans- 
formation completely decouples the charge modes of the 
IDEG from the environment. This already constitutes 
a partial solution of the problem. Moreover the results 
are generic for all values of the couplings in the basin of 
attraction of the paired-spin-liquid fixed point because, 
as we shall show, AV"c and AJc are perturbatively irrel- 
evant. 



3. Transformation to Holon Variables 

Having separated spin and charge, it is useful for many 
purposes to express the charge excitations as spinless 
fermions, which we shall call "holons". For the environ- 
ment Hamililtonian this is accomplished by rescaling the 
charge fields of the environment by the real space version 
of a Bogoliubov transformation: 
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(36) 



which also changes Kc — > "i-Kc- Then, using Eq. ( |l3| ) for 
spinless fermions, the Hamiltonian for the environmental 
charge excitations can be writen 

.93 L7.t J. , jj^-^ (37) 



where vp — Vc{4:K^ + l)/AKc and g = 27r{jf ( ^^|_^| ). 

The holons, which are created by the operator ^, are 

free fermions at the point Kc — 1/2, or 5 — 0. We 
can similarly refermionize the charge part of the pair- 
tunnelling term to obtain, when A = 1, 



^1 "^pair Ui 



^ ira 



2,c 



H.C.] 



xcos[v27r0s] cos[v 27r(/)s] 



(38) 



Thus, the pair-tunnelling term couples the holon pair 
creation operator in the environment to the joint spin 
fluctuations of the IDEG and the environment. (In this 
way, pair-tunnelling can, under the right circumstances, 
induce a spin-gap in the environment, even if there is no 
preexisting gap.) Finally, the charge density and current 
density interactions between the IDEG and the environ- 
ment (Vc and Jc) can be written simply in terms of the 
usual fermionic expressions for the charge and current 
densities, respectively. 

A similar transformation to holon variables may be 
made for the charge degrees of freedom of the IDEG. 



IV. THE PSEUDOSPIN MODEL 

The general model discussed in the previous two sec- 
tions cannot be solved exactly, although it can be studied 
using the sophisticated mean-field theory which will be 
introduced in Section V. However, the low-energy physics 
may be extracted from the solution of any model which 
has the same degrees of freedom and symmetry as the 
original model, and is controlled by the same strong- 
coupling fixed point. Here we introduce a "pseudospin" 
model which preserves the essential physics, yet it is ex- 
actly solvable Jsst . 

The essential point is that the frustration of the motion 
of holes in an antiferromagnet ]5C| ] implies that the in- 
teraction between holes in the environment is effectively 
strong, i.e. Kc and Vc are small. Thus we may ignore 
the bandwidth of pairs of holons in the environment and 
characterize them by a single renormalized excitation en- 



ergy e*. Then we introduce a pseudospin operator, r 
such that — +1/2 if there is a holon pair in the en- 
vironment in the neighborhood of R, and = ~l/2 
otherwise. (Formally, if Kc ~ 1/2, then it follows from 
Eqs. (37) and (|25|) that the pseudospin raising operator 



is given by r+ = V'lcV'Jc-) Since the pseudospins are 



discrete variables, we must put them on a lattice, where 
the lattice constant represents the distance the holon 

can diffuse in an imaginary time 1/s* . {S,p ~ \J v"^/ Kc£* .) 
Evidently, the lattice spacing is the residual effect of the 
holon bandwidth in the environment. 

The (transformed) Hamiltonian can be expressed in 
terms of the pseudospins as 



U\HpseudoUi — HiDEG 



H, 



(39) 



/ ^ JspTji COS 



'2-Kq 



R 

-E 

R 



cos[v 27r0s] 
2y/2/^^Vcdx(t>c][Tl, + l/2] 



where Hideg is the Hamiltonian of the IDEG (with 
53 = 0) defined in Eq. (^, Hs is the Hamiltonian for 
the environmental spin degrees of frcedomLwhich is the 
environmental pieces of Hs defined in Eq. ([l7| ) , Ui is de- 
fined in Eq. (p8|), and for simplicity we have ignored the 
term proportional to AJc, which we expect to be small. 
The sum is over sites in the pseudospin array, and it is 
implicit that the terms involving the continuous fields are 
integrated over a cell of size about the site R. We will 
refer to this simplified model of the dynamics of the envi- 
ronmental charge degrees of freedom as the "pseudospin" 
model. 

It is important to note that the pseudospin model 
could have been introduced at the outset to represent 
the active environment, without reference to a more de- 
tailed electronic model. In that case, Hpseudo could be 
written in terms of the original variables as 

Hpseudo — HiDEG + Hg 

Jsp[PRr+ + H.c] cos[%/2^^s] 



-^[e + 2Kpfl][r^ + l/2], 



where 



and 



^R ^ 



PR 



\s-R\<if,/2 



dxPHa 



dxp(x), 



(40) 



(41) 



(42) 



\s~R\<i^/2 



are the pair creation and charge density operators de- 
fined in Eqs. (|ll|) and (|lj) respectively, and manifestly 
[r^ -I- 1/2] is the holon pair density operator in the en- 
vironment. To see that this is equivalent to the form 
of the pseudo-spin model discussed above, we apply the 
pseudo-spin version of the unitary transformation, Ui, 



U = exp{-iV2^^r^6'c}, 



(43) 



R 



to Eq. (^). In this way, we obtain the transformed 
version of Hpseudo given in Eq. (Eol) with e* = e — 2Vc. 
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It is clear from the derivation that Hpseudo has the same 
symmetry as the starting Hamiltonian. 

In the pseudospin model, the umklapp scattering ((73) 
term of Ti.^ is unchanged by the transformation U, since 
the argument of the cosine is displaced by the trivial 
phase 47rT^, with t|, = ±1/2. Thus, in the pseudospin 
model, the canonical transformation decouples the charge 
degrees of the IDEG from the environment, even in the 
presence of a non-zero g^l 

The pseudospin model clearly captures the essential 
physics of charge fluctuations in the environment in the 
limit of small kinetic energy. In addition it is more gen- 
eral, insofar as it is also a reasonable representation of 
the spin gap centers, discussed above. Of course a con- 
tinuous distribution of centers corresponds to the case in 
which there is an environmental spin gap everywhere. 



V. EXACT RESULTS FOR THE PSEUDOSPIN 
MODEL WITH e* = AT T = 

In this section, we present an exact solution of the 
pseudospin model, Eq. (^), at a suitably chosen de- 
coupling point, in order to elucidate the mechanism by 
which a stripe coupled to a magnetic insulating environ- 
ment by pair hopping develops a gap in its spin excita- 
tion spectrum. We treat both the case in which there is a 
preexisting environmental spin gap and the case in which 
the environmental spin excitation spectrum is gapless. In 
both cases, the ground state of the solvable model is a 
fully gapped paired-spin-liquid state. However, we con- 
sider the former case to be the more physically relevant, 
as without a preexisting environmental spin gap, it is less 
likely that the model with physically reasonable values of 
the bare interactions will lie in the basin of attraction of 
the paired spin liquid fixed point. A gapped spin liquid 
is the one-dimensional version of singlet superconducting 
pairing, although it also displays enhanced charge density 
wave correlations p6| , p0{ 

A. The decoupling limit 

The close formal relation between the pseudospin 
model Hpseudo and a Kondo lattice suggests that there 
is a counterpart of the solvable limits of the one dimen- 
sional Kondo and orbital Kondo |Q arrays that we 
have analyzed previously. This is in fact the "decou- 
pling limit", discussed earlier, in which AVc — {i.e. 
Vc = nVc/'iKc), so that the unitary transformation, U, 
decouples the charge degrees of freedom of the IDEG 
from the remaining degrees of freedom. The spin part 
of the Hamiltonian remains nonlinear and, in general, it 
involves the dynamics of the pseudospins. However, a 
further great simplification occurs in the limit e* — s- 
(i.e. £ = 2Vc) at which point the pseudospin opera- 
tors, T^, commute with the transformed Hamiltonian, 
W HpseudoU , so the set of eigenvalues, rf^ — ±1/2, are 
good quantum numbers. 



In the ground state, the transformed pseudospins are 
ordered, i.e. = Tq for all R, and there is a two- fold 
degeneracy, corresponding to Tq = ±1/2. This does not 
correspond to long-range superconducting order, (which 
is forbidden in one dimension) even though the untrans- 
formed creates charge 2. After the unitary transfor- 
mation in Eq. (p9|), becomes the gauge-invariant or- 
der parameter which characterizes the composite pairing 
of the holons, and it cannot be expressed as a local func- 
tion of the original physical fields. A similar composite 
ordering was discovered for the two-channel Kondo prob- 
lem pot - Here the only symmetry that is broken in the 
ground-state is the discrete "r" symmetry, discussed in 
Appendix B. 

We show below that, so long as Jsp <C W, the array 
of pseudospins is so dense that its discreteness may be 
ignored in the ground state . Then the spin fields are 
governed by the double sine- Gordon Hamiltonian 

its — itg + rig 
Jsp f°° 
27ra 

X cos [\/27r0s (x)] . cos [V2tt (j)s {x)] . (44) 

where and Til are given in Eqs. ( p^ and (pO|), respec- 
tively. We can obtain exact solutions of the spin part of 
problem in two different limits. 

1. The case of an environment with a large spin gap 

We first consider the case in which there is a preex- 
isting spin gap in the environment and show how it is 
communicated to the IDEG. In terms of our model, this 
corresponds to the case in which Ks < 1 and |5i/ws| is 
large. Then the term proportional to gi is relevant (in 
the renormalization group sense) and even in the absence 
of coupling to the IDEG produces a spin gap in the 
environment. At energies and temperatures small com- 
pared to As, the fiuctuations of (j)s are effectively pinned, 
and cos [•\/27r0s(a;)] in Eq. ( p4[ ) may be replaced by its 
expectation value. Thus, for large environmental spin- 
gap, we can readily integrate out the environmental spin- 
degrees of freedom, leaving us which a simplified pseudo- 
spin model in which the environmental spin-degrees of 
freedom no longer appear, but in which a new effective 
coupling constant 

Jsp = Jsp < cos[V2^4>s] > ■ (45) 

replaces J^p cos['\/27r0s] in the pseudo-spin Hamiltonian, 
Eq. (^), where < cos[\/2tt^s] > is the zero temperature 
expectation value. (This expectation value can be com- 
puted exactly in the contimmm limit, < cos[\/27r0s] >~ 
As/W, from known results for the sine-Gordon field the- 
ory, as discussed below; in the strong coupling limit, 
As^W, < cos[V2^4>s] >^ 1-) 

Once this replacement is made, the analysis of this 
equation is simplified by the fact that the gi contribu- 
tion to H( is irrelevant, provided gi is not too large: On 
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the one hand, with respect to the non-interacting fixed 
point defined by Tig, the final (pair-tunnelling) term in 
Eq. ( ^4[ ) is perturbatively relevant, while the gi term is 
pertm'batively irrelevant. More to the point, the term 
proportional to Jsp is a relevant perturbation relative to 
the full sine-Gordon Hamiltonian Hq whereas if we 

reverse the logic, and we treat the gi term as a perturba- 
tion, we find that it is irrelevant. We therefore drop the 
gi term for the present with the result that Tis is reduced 
to a (solvable) sine-Gordon Hamiltonian for the field (j)s- 
As discussed below, the solution of this problem is quali- 
tatively described by the classical limit, in that (ps is thus 
pinned in the ground state, and there is a corresponding 
spin gap. 



2. The case of small, bare environmental spin gap 

When the environment does not have a large, preex- 
isting spin gap, we may omit in Eq. (Q) , and rewrite 
Hs as 

ris = At, 



with 



J sp 

Aira 



dx 



X { cos [V^cpt {x)] + cos [%/47r07 (a;)] } , (46) 

where (j)f = (0s ± (j)s)/V^- Then, in the special case in 
which the spin Hamiltonians of the stripe and the en- 
vironment are symmetric (Kg = Kg and Vg — Vg), Jis 
may be written as a sum of two independent sine-Gordon 
Hamiltonians in the variables (j)^ . The major difference 
from the case in which the environment has a spin gap is 
that Ks is replaced by 2Ks. 



B. Sine-Gordon models 

Until now, we have considered in parallel the cases in 
which the environment has and does not have a preexist- 
ing spin gap. To streamline the subsequent discussion we 
will focus solely on the more physically interesting case in 
which there is a large preexisting environmental spin gap; 
the other case can be straightforwardly analyzed along 
similar lines. So, for example, the double sine-Gordon 
model in Eq. ( |44| ) will be replaced by the ordinary sine 
Gordon model in which jTsp replaces Jsp cos [\/27r0s]. 

The solution of the resulting sine-Gordon Hamiltoni- 
ans is well known The excitations are massive soli- 
tons (which correspond to a "magnon" with spin 1 and 
charge 0) with energy spectrum given by 



where 



As 



3^ sp(^ 



(47) 



(48) 



(49) 



provided Kg < 4. In addition, so long as a < 1, there are 
breather modes [|69| , i.e. two magnon bound states, with 
spin zero and energy ~ As. In particular, as discussed 
in Eq. (|2^), spin rotation invariance implies that at low 
energies, « 1, which in the case of a large environ- 
mental spin gap implies a ~ 2/3. For a — 2/3 there are 
two breathers, one with energy As and the other with 
energy a/SAs. The spin gap As also defines a correlation 
length, ^s — I's/As, which characterizes the response of 
the spin field to external perturbations. Clearly, it is con- 
sistent to ignore the discreteness of the pseudospin array 
so long as ^s > Cp- 

There are two other classes of excitation of the spin 
degrees of freedom, both of which are non-propagating 
in the decoupling limit, but which acquire a finite (but 
large) mass when perturbations are included. The first 
involves a kink in the pseudospin order, so that, for in- 
stance, = 1/2 for < and = -1/2 for > 0. 
This induces a corresponding "half" soliton in the (f>s 
field, and so corresponds to a "spinon" with charge 
and spin 1 /2 with a creation energy, 



Aspmon ^ As 



(50) 



it is unclear at present whether 2Aspmon is greater than 
or less then As , which ultimately determines whether the 
magnon is stable or subject to decay into two spinous. 
(Classically, i.e. in the Kg —> limit, 2 Aspmon = 
\/2As > As.) The second excitation involves a flip of the 
pseudospin at one point | |70| . Again, because the spin 0s 
fields are quite rigid {i.e. ^s is large), they will hardly 
respond to such a flip, so the energy of this excitation 
can be estimated as 



S — (Jsp/TTfl) < cos(\/27r0s) > 



(51) 



(The fact that this excitation involves minimal relaxation 
of 0s can also be seen, a posteriori, from the fact that 
^< As.) 



C. Correlation Functions 

Since a continuous symmetry cannot be broken in one 
dimension, the "state" of the system is characterized by 
the correlation functions of the various possible order pa- 
rameter fields. 

In the case of noninteracting electrons, density-density 
correlation functions decay as l/cc^. Therefore, any corre- 
lation function Ci{x,x') =< Oi{x)Oi{x') > which decays 
as is "enhanced" if Ui < 2; the corresponding sus- 
ceptibility diverges as T"~^ in the limit T ^ 0. The or- 
der parameters whose correlation functions arc enhanced 
are: 2kF charge density wave 



O 



CDW 



[V'2,tV'i,t 



(52) 
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and singlet pairing 



Osp = pH^ 



(53) 



where is defined in Eq. (|ll|). At temperatures small 
compared to the spin gap, A^, the spin field is massive so 
the spin fluctuations contribute a multiplicative constant 
to these correlation functions, while all others exhibit ex- 
ponential decay. Away from half filling, there is a band of 
solitons and the exponents are given by ac dw = K* and 
asp = iZ-ft'c- Here, is the value of Kc, renormalized 
by umklapp scattering. 

For 1/2 < < 1, both SP and CDW correlations are 
enhanced, but the CDW correlations decay more slowly 
with X. However, as usual for quasi one-dimensional sys- 
tems, disorder and the coupling between stripes deter- 
mine the fate of an array of stripes. 

Even at zero temperature, the correlation function of 
the untransformed pseudospin operators decays rapidly 
with distance. However, the transformed pseudospins 
< Ut^t^'U^ > exhibit long-range order at T = 0, and 
Ising-like behavior at finite temperature, 

< U^T^nrR'U >~ (mO'cxp[-|i?-i?'|/^r(T)], (54) 

where the temperature dependent values of ^^-(T), which 
diverges as T — s- 0, and m^, which approaches 1/2, are 
estimated below. As in the case of the quantum Hall 
effect fflll. integer spin chains [Q, and various Kondo 
arrays 1^,^, so also in the present case the coherent 
state of the system is characterized by the long-range 
order of a non-local order parameter. 



VI. APPROXIMATE RESULTS FOR THE 
PSEUDO-SPIN MODEL AT T > 

Our purpose in this section is to obtain a more com- 
plete (but approximate) solution of the model at finite 
temperature and finite e* . We will also discuss, quali- 
tatively, the perturbative effects of deviations from the 
decoupling limit of the model {i.e. the effects of non-zero 
AVc). Again, for simplicity, we restrict our attention to 
the more physically interesting case in which there is a 
large preexisting environmental spin gap; the other case 
can be straightforwardly analyzed along similar lines. 
Recall that in this case, the environmental spin degrees of 
freedom can be integrated out leaving us with the pseudo- 
spin Hamiltonian, Eq. (^0[), with the effective coupling 
jTsp, defined in Eq. (|45|), replacing J^p cos[\/27r0s]. 

(It is also important to remark that the general model 
considered previously can be treated at the same level of 
approximation. The results differ little from those we ob- 
tain, here, for the pseudospin model, which substantiates 
our view that there is little physically important differ- 
ence between the two models. However, we have been 
unable to obtain analogs of the exact results discussed in 
the previous section for the general model.) 

We have shown in the previous section that the trans- 
formed pseudospins are condensed at T = 0. The impor- 
tant thermal fluctuations which destroy this order are 



the spinon excitations which produce kinks in the order 
parameter field, as discussed above. Thus, the trans- 
formed pseudospin correlation functions at low temper- 
ature are equivalent to those of a classical Ising model 
with exchange coupling, Aspinon- As a consequence, for 
sufficiently small T, the correlation length diverges as 



~ ^^J^{A spinon / P] • 



(55) 



At first blush, Eq. (|55|) might be expected to apply 
so long as r <C Aspmon, but in fact it only holds so 
long as T ^ S; this is because at temperatures of or- 
der 6, the large density of thermally excited single pseu- 
dospin flips (which, by themselves, directly affect only the 
magnitude, but not the range of the pseudospin order) 
leads to a large renormalization of the spinon creation 
energy; Eq. (|55| ) remains valid, but with a temperature 
dependent renormalized spinon creation energy replacing 



^spinon 



(and latice constant, ^p). 



We obtain an estimate of this renormalization using 
the technique of Coleman, Georges, and Tsvehk Ba- 
sically, this amounts to making a mean-field like decom- 
position of the non-linear term (i.e. the term propor- 
tional to J'sp) in Hpseudo, so that in computing the ther- 
modynamic properties of 4>s , we replace the transformed 
pseudospin operators by their thermal expectation 
value, ttIt- U^rf^U >, and conversely in computing 



the pseudospin properties, we treat < cos[v27r0s] > as 
a pseudo magnetic field. As with all mean-field theo- 
ries in one dimension, this approximation has the fault 
that it produces spurious long-range order at finite tem- 
perature, where < U^rf^U > and < cos[^/2TT(j)s] > are 
actually equal to zero. However, we shall see that the 
mean-field theory is exact in the limit e* and T tend to 
zero, and thus its results are reliable at low temperatures 
when it is used to estimate local quantities such as A^, 
Aspinon and m^. In other words, it is correct for inter- 
mediate scale fluctuations. (For example, nir should be 
deflned in terms of the asymptotic form of the composite 
order parameter correlation function in Eq. ( |5^ ) , and the 
mean-field theory should be viewed as a way of estimat- 
ing it as the "local" expectation value of an operator.) 

In the mean-field approximation the self-consistent 
equations for the temperature-dependent gaps, As{T) 
and S{T) are: 



6{T) = (2Jsp/7r) < cos[V2^0, 
As(r) = As[2<r- >]", 
< C/V^f/ >= {6/4Eb)taiih[f3Eb{T)], 
Eb{T) = vW + W2F, 



>, 



(56) 
(57) 
(58) 
(59) 



where As and 5 are, respectively, the values of As and 
5 at T = and e* = 0, as given in Eqs. ( ^ ) and 
( ^T| ) above. Finally, < cos[\/27r0s] > should be com- 
puted at finite temperatures using known results from 
the thermal Bethe ansatz [|7^ for the sine-Gordon model. 
These results are quite complicated, but fortunately the 
information we need is fairly minimal; specifically, that 
< cos[\/2TT(f)s\ > is a monotonically decreasing function 
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of temperature, with the scale for the temperature de- 
pendence set by the zero-temperature gap. Among other 
things, this imphes that so long as As(T) ^ T, we can 
use the zero-temperature result 



< cos 



> 



(2Q-1) 



(60) 



for the sine-Gordon part of the calculation. It is clear 
from these equations that, for T <C £'b(0) and T <C 
As(0), all gap parameters are well approximated by their 
zero temperature values. Conversely, the gaps begin to 
decrease when for T ~ -^^(O) if Eb{0) < As(0), or when 
T - As{0) if A^(0) < Eb{0). We can, in general, define 
a characteristic crossover temperature, Tpair, to be that 
temperature at which As(T) begins to drop significantly 
from its zero-temperature value. In some cases, this is the 
only obvious crossover temperature in the problem. How- 
ever, we will see that under some circumstances, it is still 
true that As (T) ^ T for a substantial range of tempera- 
tures above Tpair', in these cases there is a second, para- 
metrically larger crossover temperature, T^^^^ ^ Tpair, 
at which the spin gap gets to be comparable to T. For 
temperatures above T^^^^, all effects of pairing and co- 
herence are negligible. 



We can now proceed to analyze the solution of these 
equations as a function of temperature and e*. The re- 
sults (for the important case mandated by spin-rotation 
invariance in which a = 2/3) can be sumarized as follows: 
As(0) is largest for e* = 0, and falls slowly, roughly as 
£^^, with increasing e*, but only vanishes {i.e. pairhop- 
ping becomes irrelevant) when e* [JspY /gi- Tpair is 
much smaller than A^ (0) for small e* , but increases with 
increasing e* , reaching a maximum for e* Jsp, at which 
point all energy scales are comparable; Tpair ^ As(0) ~ 
Jsp. Meanwhile, T'^^^^ is of order J^p and roughly inde- 
pendent of e* for £* small compared to jTsp, and becomes 

> Jsv These results 



indistinguishable from Tpair for ^* 
are shown schematically in Fig. H. 



In the following, we derive these results, focussing se- 
quentially on four distinct regimes of behavior as a func- 
tion of £*; in the subsection headings, the ranges are ex- 
pressed with numerical exponents for the important case 
a = 2/3, as well as algebraically for general a. 



log [Energy/e?.,] 




log [A./e'.p] - 


J. 


A.{0) 


- 












log [8/e'sp] - 




1 ' 



log [S/c»,p] 



log [e*/cy 



FIG. 2. Energy scales from the solution of the pseudospin 
model as a function of e*: & and As are, respectively, the 
coherence scale and the spin gap derived from the exact solu- 
tion of the model for e* = and given in Eqs. (48) and (51), 
As(0) is the zero temperature value of the spin gap, Tpair is 
the temperature scale at which Ai,(r) begins to fall signifi- 
cantly relative to its zero temperature value, and Tp^i^ is the 
temperature at which As(r) ~ T. 



A. For the case e* < Jsv{JsvlWfl^ , 
i.e. when e* ^ 5 

In this regime, the results are qualitatively the same 
as for £* = 0. (Note_that for £* = T = 0, the self- 
consistent equations 



60) are exact.) There is little 



temperature dependence of any of the gap parameters in 
the low temperature regime, T <^ Tpair, where 

Tpair ^ ^- 



(61) 



Clearly, substantial suppression of As(T) due to pseu- 
dospin fluctuations begins to occur when T ~ Tpair', as a 
consequence, Tpatr/^siO) ~ piEf)As{0) < 1/2. 

There follows an intermediate temperature regime. 



tpajr < T < T^air^ whcrC 



T' 

pair 



Ss[A/l 



l2(l-Q)/(2-a). 



(62) 



in this regime, even though As{T) is strongly suppressed, 
it is still true that As{T) 3> T, so we can approximate 
< cos[\/2Tr (j)s] > by its zero temperature value Eq.(|60|), 
with the consequence that 



AsiT)^A,[f36] 



V2(l- 



and 



(63) 



(64) 
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However, while significant spin pairing still survives 
in this temperature range, the entropy of the pseu- 
dospins is recovered, and hence the specific heat, Cy ~ 
[(5(0)/r]i/(i-"), is large. 

'^pair is the temperature at which T — As(r), where 
As(T) is given by Eq. (p3|). For temperatures T ^ 
T^air^ there is no coherence, no apparent gap in any of 
the degrees of freedom, and the problem can be treated 
using a high temperature expansion. 

We can summarize the heirarchy of scales in this case 

as 

A, ~ A,(0) » r;„„, » Tparr ~ m ^ ^ » £*. (65) 

Specifically, for the a = 2/3 case. As ^ ^Jsp , ^a^r ^ 
\J 3p ; and S ^ Jsp ■ 

B. For the case Jsp[Jsp/W]^^^ < e* < Jsp 

I.e. when J < £* < 5^ [A./^sl^f^-^^/'^-"' 



C. For the case Jsp < e* < W', 

I.e. when ^[As/^]^^!-")/'^-") < e* < W 

Whenever 5[A < e*, it follows that 
As(0) <C e*. As a consequence, the temperature depen- 
dence of the various gaps is set by 

Tpa^r - A,(0) (70) 

where A3(0) and S{0) are given by Eqs. ( ^8|) and (|6^), 
above; moreover, there is no longer a distinct tempera- 
ture scale Tp^j^. 

The heirarchy of scales in this case can be summarized 

as 

A,> A,(0) - Tpa^r > ^ > S{0), 

e*» A,(0). (71) 

In this regime, both As(0) and, correspondingly, Tpair are 
decreasing functions of e*. To be specific, f or the case of 
a = 2/3, Tpa^r ~ J^p/e* and ,5(0) - Tp^.^yFTW^. 



It is easy to see from Eqs. (|8|) and ( |5^ ) that larger val- 
ues of £* supress the thermal disordering of the pseudo- 
spins, and hence removes the anomalous renormalization 
of As(r) at low temperatures. At T = 0, and so long as 
e* » S, 

5(0) =^[^/e*]i/2(i-«) (66) 

and 

A,(0) = A3[^/£*](2"-i)/2(i-"). (67) 
If at the same time, e* < Tpair' then As(0) > e*, so 

Tpair ~ £*■ (68) 

For T <C Tpair, there is little temperature dependence 
of the gaps, whereas, for T ^ Tpair, £* falls out of the 
problem so As{T), S{T), and T^^j,. are given by Eqs. 
(g), (H), and @, as before. 

The remarkable property of this range of parameters 
is that, as e* increases, the spin gap at T = decreases 
rapidly (as expected) but the pairing temperature, Tpair 
actually increases. In other words, in order to obtain a 
high temperature scale for pairing, the charge transfer 
energy, e* , should be somewhat above the Fermi energy! 

We can summarize the heirarchy of scales in this case 

as 

A, » A,(0) » r;„„, » Tpa.r ~ £* » ^ » <5(0). (69) 

One remarkable feature of this result, which relies on 
the particular value a — 2/3, is that in this regime 
A,(0) ~ [Jsp]ye*, Tpa^r ~ s* , and T;„^ ~ Jsp are all in- 
dependent of the bandwidth. Note that at the upper end 
of this range, A,(0) - Tp„, ~ T^^„ ~ £* ~ Jsp\ This 
same conclusion follows from evaluating the expressions 
in the next subsection at the lower limit of the stated 
range. 



D. e* ~ W: Renormalized interactions 

In the limit of large e*, the dynamical nature of the 
collective mode is unimportant; it could have been in- 
tegrated out to obtain new effective interactions in the 
IDEG, with retardation and spatial non-locality limited 
by the size of e* . Moreover, since in this limit, holon 
pairs in the environment exist only as dilute, virtual ex- 
citations, it is sufficient to compute these interactions 
perturbatively in powers of Jspjs* ■ To second order in 
Jsp, the Hamiltonian is of the same form as Hideg in 
Eq. (^), but with a renormalized chemical potential and 
interactions: 

g*i^ gi- 5g. (72) 
K: = K{gl) (73) 

vl ^Vs+ <5.g/27r, (74) 

where 6g = (Jsp)^/4e*. 

When gi is small, g\ < Q and the pair fluctuations 
produce a net attractive interaction in the spin degrees 
of freedom, which leads to a spin gap of magnitude ||7^ 

A, = A^2^{vs/a) exp[-l/A] (75) 

where A = Ps\gi\/0' and Ps = a/irvs- It is also clear 
that there is a corresponding crossover temperature, 
Tpair ~ As/2 <C £*, abovc which the spin gap vanishes 
and the spin excitations are well described as linearly 
dispersing collective modes with velocity v* . Again, the 
charge modes are completely unaffected by the pairing 
physics, and so continue to be described as linearly dis- 
persing modes with velocity Vc. Hence the Drude weight 
(or, equivalently for the IDEG, the zero temperature su- 
perfluid phase stiffness ) is unrenormalized. 
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This analysis is strictly correct only if e* > W, be- 
cause it did not take account of retardation, which im- 
plies that the induced interaction Sgi, vanishes for en- 
ergy exchange much greater than e*. However, for the 
physically more interesting case, W ^ e* Jsp, the 
effect of retardation can be studied using an energy shell 
renormalization group scheme, as in the electron-phonon 
problem [po[ . This improved treatment produces results 
that are similar in spirit to those described above, ex- 
cept that, for energies smaller than e* (when there is no 
longer a distinction between the retarded and instanta- 
neous pieces of the interaction) , the effective interaction 
has a renormalization, Sgi, which is a complicated, but 
calculable (20), function of gi, (Jsp)2/4e*, and e*/W. In 
all cases, there is a critical value of the charge transfer 
energy Sc ^ [JspY I ^gi, such that for larger e* > Sc, the 
renormalized value of gi is positive at low energies, and 
there is no spin gap, whereas for e* < ec, <?i is negative 
and a spin gap opens up at zero temperature. This an- 
swers the question of how "active" the environment must 
be. 



E. Effects of "irrelevant" interactions 



Tpair ~ where strong pairing sets in. In higher dimen- 
sional systems this large entropy is presumably respon- 
sible for heavy- fermion behavior in the model Q; in the 
present context it leads to the small ratio of Tpair /^{O). 
When e* > S, the majority of the entropy associated 
with the pseudospins is lost at temperatures greater than 
Tpair- As a consequence, thermal disordering effects are 
relatively less severe, and Tpair/^{0) ~ 1/2 is rapidly 
restored. 



VII. THE BEHAVIOR OF THE CHARGE 
DEGREES OF FREEDOM: 

We have seen that, in the pseudospin model, the 
canonical tranformation decouples the charge degrees 
of the IDEG from the environment, and their fluctua- 
tions are described by the quadratic Hamiltonian, TiJ?. 
This Hamiltonian describes a fluctuating superconduc- 
tor, with phase 9c, or in dual language, a fluctuating 
charge density wave, with phase (pc- Evidently, proxim- 
ity to commensurability or the existence of an external 
potential can substantially modify the physics. 



We now consider the effects of various interactions that 
we set equal to zero in the decoupling limit. Because the 
spectrum of the pseudospin model has a gap at the solv- 
able point, all of the omitted terms are formally irrele- 
vant in the renormalization group sense. Of course this 
does not give us license to completely ignore these terms; 
they can have large quantitative, and at times qualita- 
tive effects on the physics of interest, even if they do not 
affect the character of the true asymptotic behavior of 
the system. 

Let us consider the effects of non-zero AV and, e* 
on the nature of the excitations of the system at zero 
temperature. When these couplings are small, their most 
important qualitative effect is to induce dynamics for the 
pseudospins. In the presence of these terms, the effective 
Hamiltonian for the pseudospins, obtained by integrating 
out the electronic degrees of freedom , is qualitatively 
similar to (but not precisely equal to) the spin 1/2 Ising 
model in a transverse magnetic field, 

i/^//^^P/2)r-+eV^j 

R 

- J2 [K{R - R')r^Rr^B' + K{R - i?>^r^,] (76) 

R>R' 

in which if (i?-i?') - ^VA,. £ind k{R-R') - {AV f/A, 
and both have range of order ^s- As is well known, 
a transverse field induces dynamics (propagation of the 
kinks) in the spin 1/2 Ising model. 

As we have seen, the other effect of e* is to supress 
thermal fluctuations of the pseudospins. At high tem- 
peratures, there is an entropy density S — (a/^p) In 2 as- 
sociated with the discrete symmetry of the pseudospins. 
For e* = 0, this entropy is lost at about the temperature 



A. The role of Umklapp scattering 

The charge fields of the IDEG are governed by the 
Hamiltonian: 



He — + HI, 



(77) 



where Hq and iff are given in Eqs.(p^ and (^). Now 
the c-number (4^^ — G)x may absorbed into the phase 
0c, without changing the commutation relations and the 
quadratic part of H^ in Eq. (|2^) may be diagonalized 

1/2 

by the canonical transformation (pc ipcKc , He — *■ 

1 /2 

Ilc/Kc ■ The net result is that the charge degrees of 
freedom are described by a sine-Gordon model with a 
chemical potential fi* given by 



i* = 



Vcj Akp - G) 
iKc 



(78) 



For the strongly incommensurate case, in which fi* is 
large, we can ignore the umklapp scattering term (pro- 
portional to 33 ); in this case the charge excitations are 
gapless collective modes with a sound-like dispersion and 
a velocity, Vc, that is unrenormalized by the interac- 
tions with the environment. Correspondingly, the Drude 
weight, or superfiuid phase stiffness (which cannot be dis- 
tinguished in one dimension in the absence of disorder) 
are also unrenormalized. 

In the nearly commensurate case, which character- 
izes the doped-insulator region, the analysis of the corre- 
sponding sine-Gordon theory is the same as for the spin 
degrees of freedom. In particular, for Kc < 1, which is 
always satisfied for repulsive interactions, the "particles" 
in the theory are massive solitons with charge e and spin 
0. It follows at once that the system undergoes an insula- 
tor to metal transition at = Ac, where the chemical 
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potential moves out of the gap, and that there is a finite 
density of sohtons: 



nsoi = 



A2 



(79) 



with given in Eq. ([Tq). For smaU Usoi, the Drude 
weight of the stripe is proportional to Hgoi. This ar- 
gument is similar to the analysis of the commensurate- 
incommensurate transition by Pokrovsky and Talapov 
ifzsl , except that they considered a two-dimensional clas- 
sical problem, equivalent to the quantum sine-Gordon 
problem in imaginary time. 

For quarter-filled stripes |Q, Akp = 2kp — G/2, so 
the charge density on the stripe and in the environment 
may jointly lock to the lattice. This commensurability ef- 
fect competes with superconductivity but, if the coupling 
constant is not too large, it may not develop beyond the 
logarithmic temperature dependence that characterizes 
the early stages of renormalization j?^]. We are investi- 
gating this behavior as a potential source of the special 
stability of quarter-filled stripes for doping x < 1/8 in 
the La2Cu04 family [ p|j7^ ], and the logarithmic tem- 
perature dependence of the resistivity observed [ [79|j4l| 
when the onset of superconductivity is suppressed. 



B. External Periodic Potential 

Here it is assumed that there is an external potential 
with a wave vector q which is close to 2kp. Then the 
Hamiltonian must be supplemented by a contribution 

/oo 
dxiV-I ,>2,.e^<2'=-"~''"^ + B.C.] , (80) 
-°° 

which may be written in the boson representation (|l^ ) 
as: 



2u , , 
He — — / dx cos 
na 



27r0c + (q — 2kp)x 



cos 









(81) 



It is straightforward to show that, when the pseudospin 
representation is introduced for the charge degrees of 
freedom of the environment, then He is not changed by 
the unitary transformation defined by [/ in Eq. (28): i.e. 
W HeU = He. Moreover, it is clear from the spin Hamil- 
tonian ( |44| ) that cos[\/2TT(j)s{x)] has a finite expectation 
value so that it may be replaced by a constant in He 
to obtain the asymptotic behavior of the charge degrees 
of freedom. Umklapp scattering may be ignored if it is 
an irrelevant variable or if 4:kp is sufficiently far from a 
reciprocal lattice vector. However, the effect of the pe- 
riodic potential is similar to that of umklapp scattering. 
The main differences are that the solitons are massive 
when Ke < 4, (as opposed to Ke < 1 for umklapp scat- 
tering) and that fi* = Vc{2kp — q)/Kc, which modifies 
the condition for the metal-insulator transition. 

The physical argument for including such a potential is 
as follows: In the ordered state of Lai.6-3;Ndo.4Sr^Cu04, 



the holes on a given stripe move in an effective potential 
produced by the stripes in a neighboring Cu02 plane. 
Since stripes in adjacent planes are perpendicular to each 
other, the wave vector of the charge contribution to the 
effective potential is given by g = 2e in units of 27r/a, 
where a is the lattice spacing In the same units, 

2kp — ns/2, where is the concentration of doped- 
holes on a given stripe. The present experimental evi- 
dence 1 41 4^ is consistent with e = 1/8 and — 1/2 
and hence q = 2kp for dopant concentration x = 1/8. 
This is the hole concentration near which the supercon- 
ducting Tc is suppressed in the stripe-ordered material 
Lai.6-a;Ndo.4Sra;Cu04 Q and in La2-a;Baa;Cu04, for 
which there is indirect evidence of stripe order | pl[ | . 

An array of stripes will undergo a transition to a super- 
conducting state at a temperature which is determined 
by the onset of phase coherence and is proportional to 
the superfluid phase stiffness j|] which, in tm'n, is pro- 
portional to Hsol- 

In Sec. HI we considered the case in which the en- 
vironmental spin gap arose because the backscattering 
term proportional to gi was relevant. For a half-filled 
band with ^3 also relevant, there is a broken symmetry 
ground state with period 2a, which produces an exter- 
nal potential on the stripe, with a wave vector equal to 
Akp when Us = 1/2. Such a potential is commensurate 
with the umklapp term (73 , so the coupling between these 
terms must be taken into account. This is an example 
in which spin gaps with and without a broken symmetry 
may lead to different consequences. The physical case 
has no broken symmetry. 



VIII. SPIN GAP CENTER 

Another model of some physical interest has a spin gap 
at one specific location as, for example, at an isolated an- 
tiferromagnetic region in a metal. This is an example of 
a dynamical impurity problem, in which the conduction 
electrons couple to a center with internal degrees of free- 
dom. It is well known that an angular momentum analy- 
sis produces a one-dimensional Hamiltonian involving the 
radial motion of incoming and outgoing fermions on the 
half line r > 0, where r is the distance from the pairing 
center |6^. Also, it is possible to extend the space to all 
values of r by transforming incoming fermions for r > 
to incoming fermions at position — r. Then the problem 
is formally equivalent to a one-dimensional electron gas 
in which only the right-going fermions interact with the 
pairing center. In the absence of left-going fermions, the 
operator P\ introduced in Eq. (pT|), cannot be defined 
and only the r^-pairing term I84I 



(82) 



couples to the pairing center. Triplet pairing terms are 
omitted because the exclusion principle requires them to 
be of the form ip\ | which is less relevant than j^. 
(The derivative in the triplet operator leads to an extra 
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power oil/x^ in the correlation function.) Thus a pairing 
center naturally produces singlet pairing. 

We consider the case in which the center has a large 
spin gap, so the pseudospin variable (representing charge 
transfer to the center) is the only internal degree of free- 
dom of the center that we retain, explicitly. Thus the 
Haniiltonian is 

^center — tl\ 

DEG + (83) 

where Hideg is given in Eq. (|2|), although in the case in 
which the metallic degrees of freedom represent a higher 
dimensional Fermi liquid, one must set the interactions 
{go) to zero. The bosonized form of is 



V 



y27r 

7ra 



(84) 



Here $i,c(0) = [$i,t(0) + $i,x(0)]/\/2. In this form the 
model is equivalent to a single-channel Kondo problem 
pO| , and it may be solved by making a unitary transfor- 
mation H center ^ U^HeenterU with 



U = ea;p[-a$c(0)r^] 



(85) 



and choosing A = V2 — 1, for the special point V ~ 
•\/27rAuc- Then H center becomes 



W HcenterU— HiDEG + 



H.c] 



(86) 



This the Hamiltonian may be "refermionized" by writing 
the pseudospin operator in the form t+ = -qd where rj is 
an anticommuting c-number and d is a fermion annihila- 
tion operator, and inverting the boson representation of 
fermion fields: 



■01 = 



(87) 



When written in terms of these variables, the right-going 
part of the Hamiltonian becomes 



ivc / dx[ipldxipc 



:(dV4(0) + i/.C.: 



which is precisely the Toulouse limit from which all of the 
well-known behavior of the single channel Kondo problem 
may be derived ]60| . This argument strongly suggests 
that arrays of pairing centers in two and three dimensions 
behave like Kondo lattices, and that they should show 
heavy- Fermion behavior 

Of course a single pairing center in a purely one- 
dimensional model should also exhibit this single-channel 
Kondo behavior. This would not happen if we replaced 
the pseudospin array in Eq. ( pO| ) by a single center, be- 
cause we would have omitted a possible 77-pairing interac- 
tion, of the form j/,jT^[P^ 1 -t- Pri.2] in that Hamiltonian. 



While momentum conservation indeed makes this term 
unimportant for the extended array, a spin-gap center, 
by its very nature, breaks translational symmetry and 
hence permits finite momentum transfer scattering pro- 
cesses. Including these terms, the total pair coupling at 
a single spin-gap center in Eq. (|l^) may be written 



Upair — '(i^sp-P^(^) 

+ j;,[i^„.i(i?)e'^'="^ 



ri.2{R)e' 



(89) 

H.C. 



If we consider a single center at i? = 0, and consider 
the case jTi, = 0, the left-going fermions at position x 
may be transformed to right-going fermions at position 
—X, without changing the Kondo coupling. Thus the 
subscripts 1,2 become "flavor" labels and we have a two- 
channel Kondo problem. However, in this language, the 
Jri term breaks the "channel degeneracy" , and is pertu- 
batively relevant, so it produces a single-channel Kondo 
problem. On the other hand the oscillating factors in Eq. 
( p9| ) make the Sri perturbatively irrelevant for the array 
and moreover, since the mismatch of momenta between 
the IDEG and the antiferromagnet imply that is small 
compared to J^sp , the neglect of 7y-pairing interactions for 
the extended system is justified. This is analagous to the 
behavior found previously for Kondo systems p9| , where 
the anisotropic single-channel Kondo array, behaves as 
if it were a two-channel Kondo array, even though the 
single impurity version of the model exhibits ordinary 
Kondo behavior. 



IX. DISCUSSION 
A. Summary of results 

We have studied a model of a IDEG in an active envi- 
ronment, focussing in particular, on the case in which the 
environment possesses both a charge gap and a spin gap, 
and the energy difference between a singlet pair of holes 
in the IDEG and the environment, e, is small in com- 
parison to the bandwidth. We have discovered a new 
mechanism for producing strong superconducting fluc- 
tuations on a high temperature scale, in which a spin 
gap is induced in the regions between the stripes by spa- 
tial confinement, and transferred to the IDEG by pair 
tunnelling. A striking feature of this mechanism of su- 
perconductivity, which may be described as a spin-gap 
proximity effect, is that the pairing (i.e. the spin gap) is 
a property of the insulating state itself, and it is simply 
imprinted on the mobile holes through their virtual ex- 
cursions into the insulating regions. We have found that 
this phenomenon is robust and, in particular, it survives 
the presence of strongly repulsive forward scattering in- 
teractions, i.e. Coulomb repulsion between electrons. 

We have demonstrated that the physics of this prob- 
lem is captured by a simple pseudospin model, for which 
exact and well-controlled approximate results can be ob- 
tained. This model includes the most important interac- 
tions: the renormalized pair-tunnelling matrix element, 
Jsp (defined in Eq. (|45|)), the renormalized energy cost. 
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£*, required to move a singlet pair of holes from the 
IDEG to the environment, the bandwidth of the IDEG, 
W ~ Ep, (which is assumed to be large compared to 
other energies), and the exponent, a, which character- 
izes the spin correlations of the IDEG. We have used 
renormalization group arguments to show that a « 2/3 
for repulsive, spin-rotationally invariant interactions, and 
we shall use this value of a in discussing our results. 

We have found that, generically, this model produces 
singlet pairing (spin gap behavior) at a high tempera- 
ture, Tpa.r- in the limit e* 0, Tpair ~ JspiJsp/W)^^^ , 
while for e* > JspiJsp/W^/^ , Tpair is the smaller of 
e* and As(0) ~ J'^pjs* . Remarkably, this means that 
for small e*, Tpair is an increasing function of e*, which 



reaches a maximum value of T, 



pair 



Jsr, when e* ^ J, 



sp • 



Below Tpair, singlet superconducting and CDW suscep- 
tibilities diverge as T ^ 0, with the stronger divergence 
typically associated with the CDW. Moreover, this high 
pairing scale is not accompanied by any significant re- 
duction of the zero temperature superfluid phase stiff- 
ness (Drude weight), i.e. there is no strong mass renor- 
malization. We have also identified a zero-temperature 
spin gap energy, As(0), which plays the role of the su- 
perconducting gap, Aq. In the small e* limit the ra- 
tio Tpair/^s{0) ~ ^s{0)/W « 1/2, while for large e* , 
Tpair /As{0) « 1/2, as in BCS theory. (The evolution of 
these energy scales as a function of e* is shown in Fig. ||, 
and discussed in Sec. VI.) The ground state of this model 
has a broken, discrete Z(2) symmetry, unrelated to any 
of the usual space-time symmetries of the problem, and 
a corresponding non-local order parameter which devel- 
ops a non-zero expectation value in the ground state, and 
has an exponentially long correlation length at low tem- 
peratures. (See discussion of "r" symmetry in Appendix 
B.) 



B. Interactions Between Stripes and Possible 
Ordered Phases 



To extend our results to situations in which there is 
a true phase transition, we must consider the proper- 
ties of an array of one-dimensional systems (stripes) . To 
avoid misunderstanding, we emphasise that, for purposes 
of the present discussion, "CDW" refers to charge or- 
dering along the stripe direction, whereas "stripe order" 
implies charge ordering in the direction perpendicular to 
the stripes, i.e. ordering of the stripe positions and ori- 
entations. Of course, both types of order are a form of 
generalized charge density wave. 

The ultimate nature of the long-range order de- 
pends, among other things, on the coupling between 
stripes, which is profoundly influenced by the intervening 
antiferromagnetically-correlated regions and, in particu- 
lar, by the frustration of hole motion in the antiferro- 
magnet, which was the driving force for the formation 
of the stripes themselves. Thus, this coupling should be 
smaller than the characteristic energies of the electronic 
correlations along the stripe, considered in this paper. 



With this in mind, the onset of superconductivity in a 
dilute stripe array can be studied by introducing weak in- 
teractions between well-separated stripes. Single- particle 
tunnelling between stripes is an irrelevant perturbation 
psf , because of the existence of a spin gap, so we do not 
expect a crossover to higher-dimensional Fermi liquid be- 
havior in this limit. Then the nature of the long-range 
order is determined by pair tunnelling and the Coulomb 
coupling between stripes. 

Effects of Disorder: There are two distinct types of 
disorder which have very different effects on the physics 
of an array of stripes. The first is a degree of random- 
ness in the couplings between stripes, which may be pro- 
duced by impurities (as in e.g. organic conductors) or 
by quantum or thermal fluctuations in the stripe con- 
figuration. For a "self-organized" quasi one-dimensional 
system, such as a charged stripe array, the latter source 
of disorder is likely to be the more important. Disorder 
of this type favors superconductivity (which is a fc = 
order) since it strongly frustrates the short-wavelength 
CDW order associated with the Akp or 2kF instabilities 
of the IDEG. This is especially so when the stripes are 
strongly fluctuating. In the simplest situation, the dy- 
namics of the stripes is slow compared to the Josephson 
plasma frequency, as for example in La2-a;Sra;Cu04, and 
the disorder is essentially static. On the other hand, if 
the CDW and superconducting fluctuations are on sim- 
ilar time scales, new physics may emerge; an interesting 
possibility is that there exists a novel quantum critical 
point which controls the physics in some region of tem- 
peratures and dopant concentration |]8^ . 

The second type of disorder affects the coherence of 
electronic motion along a single stripe. For a single 
stripe, the back scattering of holes from an impurity is al- 
ways pertubatively relevant for the range of interactions 
considered here, because CDW correlations are enhanced 
p5[ . However, the localization can be superseded by 
sufficiently strong Josephson coupling (pair-tunnelling) 
between stripes, and there will be an insulator to su- 
perconductor transition as the concentration of stripes 
grows or the Josephson coupling between stripes is, in 
any other way, increased, with fixed disorder. This is 
in agreement with the evolution of the ground state ob- 
served in La2-a;Sra;Cu04 as a function of doping [T^ , or 
applied magnetic field [fzoj . 

Symmetry of the order parameter: If stripe or- 
der breaks the four-fold rotational symmetry of the 
crystal, the superconducting order will have [p6|, 87 
strongly mixed extended-s and dx2_y2 symmetry! This 
will happen in a stripe-ordered phase, such as in 
Lai.6-2;Ndo.4Sr2,Cu04, or in a possible "stripe nematic" 
phase, in which the stripe positional order is destroyed by 
quantum or thermal melting or quenched disorder, but 
the stripe orientational order is preserved. (Such phases 
also would be characterized by large induced asymme- 
tries in the electronic response in the a6-plane. Below we 
discuss some preliminary evidence for a transition to a 
stripe nematic phase in overdoped YBa2Cu307_5.) 

On the other hand, when the stripes are disordered 
at long length scales, the thermodynamic distinction be- 
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tween s-wave and d-wave superconducting order is well 
defined; however, even here, if there is substantial orien- 
tational order to the stripe fluctuations at intermediate 
length scales, the interplay between the two types of su- 
perconducting order is likely to be more complicated and 
more subtle than in conventional, homogeneous materi- 
als. For example, one can imagine that, even in a phase 
which is globally c?-wave, substantial mixtures of s and 
d-wave order could occur over mesoscopic scales near sur- 
faces or twin boundaries. 

Superconducting fluctuations: A necessary corollary of 
the stripe model is that, in lightly doped materials, the 
temperature scale, Tpair, at which pairing occurs (on a 
single stripe) is parametrically larger than the supercon- 
ducting transition temperature, Tc, which is governed by 
the Josephson coupling between stripes. Moreover, since 
the pairing force derives from the local antiferromagnetic 
correlations in the regions between stripes, both Tpair 
and Tc must be less than the temperature scale, Tj^p, 
below which local antiferromagnetic correlations develop. 
A sequence of crossovers is, indeed, observed experimen- 
tally in underdoped high temperature su per conductors, 
and they have tentatively been identified ||89[| with these 
two phenomena; see Fig. above, and the discussion 
below. 



C. Phase diagram of the high temperature 
superconductors 

The schematic phase diagram shown in Fig. ^ shows 
the global framework in which our model is related to 
the properties of the high temperature superconductors. 
The axes in this figure are temperature, T, and doping 
concentration, x: hatched lines indicate the most impor- 
tant crossover temperatures, and the solid lines repre- 
sent phase transitions to the antiferromagnetically or- 
dered state at very small x, and to the superconduct- 
ing state at larger x. (In general, there are additional 
phase transitions and possibly other crossovers, but here 
we wish to focus only on the central physical issues.) 

The upper crossover temperature T^* characterizes the 
aggregation of charge (holes) into stripes; as we have 
shown elsewhere, the driving force for this crossover 
is frustrated phase separation |p^-p^. Above Ti the 
holes are more or less uniformly distributed, and ran- 
domly disrupt antiferromagnetic correlations, while be- 
low T^ , the self-organized stripe array allows local an- 
tiferromagnetic correlations to develop in the hole-free 
regions of the sample. At short distances, low energy 
spin fluctuations should come from regions with the char- 
acter of odd- leg ladders, and be like those of the one- 
dimensional Heisenberg model |^2|, and, indeed, there is 
experimental evidence poj indicating that this is the case 
in La2-a;Sr2;Cu04. As a; — > 0, T^ approaches the tem- 
perature T^ at which local antiferromagnetic correlations 
develop in the undoped systems [Q. Between Tj* and 
the superconducting transition temperature Tc, there is 
a large range of temperatures in which there are signif- 
icant stripe correlations, but coherence between stripes 



can be largely ignored; this is the region of tempera- 
tures addressed by the calculations in this paper. As the 
concentration of holes increases, the separation between 
stripes eventually becomes comparable to their width, at 
which point all information concerning the Mott insulat- 
ing state is lost; for this reason, we have shown Tj* 
at a dopant concentration, Xmax- 

We identify the lower crossover temperature Tj* with 
Tpair , the temperature at which pairing (spin gap) behav- 
ior emerges within a stripe. This is also the temperature 
below which significant local, quasi one-dimensional su- 
perconducting fluctuations become significant. For local 
probes of the spin and quasiparticle response functions, 
the system should appear all but superconducting below 
this temperature. Since Tpair is more or less a property 
of a single stripe, we have shown it as a relatively insensi- 
tive function of x, until it is cut off by Tj* at larger dopant 
concentrations. From this figure, it is clear that Tpair is 
substantially greater than Tc throughout the underdoped 
regime, and possibly even at optimal doping, and only ap- 
proaches closely to Tc in the overdoped regime. Thus, in 
underdoped materials, Tc is determined by the superfluid 
phase stiffness, and hence by the Josephson coupling be- 
tween stripes, rather than by the pairing scale. This is 
consistent with our previous analysis [||. 

It should be noted that a phase diagram of the same 
form as that shown in Fig. ^ has been considered, previ- 
ously, on purely phenomenological grounds , with the 
crossover temperatures determined as follows: 

1) The upper crossover occurred at a characteristic tem- 
perature deduced by Batlogg and coworkers from an 
analysis of susceptibility and transport properties, and 
by Loram and coworkers ||9^ from an analysis of ther- 
modynamic data. We feel that all of these phenomena 
are broadly consistent with our identification of Ti with 
the emergence of stripe and local antiferromagnetic or- 
der. (It appears that a pseudogap appears in the c-axis 
optical conductivity at this temperature. Much of 
the c-axis optical oscillator strength will be shifted to 
energies higher than + e*/2 as the stripe correlations 
emerge below Tj*.) If we accept this identification, then 
for moderate doping concentrations, a typical value is 
Ti ^ 30QK, although it depends somewhat on the par- 
ticular material, and rather more strongly on the dopant 
concentration. Indeed, stripe correlations have been seen 
neutron scattering experiments all the way up to 300K, 
although the scattering cross section decreases continu- 
ously, making it difficult to identify them unambiguously 
at high temperatures pO| . 

2) The lower crossover was identified by Batlogg and 
Emery as the characteristic "pseudogap" tempera- 
ture, deduced from the temperature dependence of the 
Cu NMR 1/TiT, which correlates well with the emer- 
gence of superconducting gap structure in ARPES ex- 
periments [^,^, and a narrowing of the "Drude-like" 
peak in the optical conductivity in the afo-plane ||96|| . If 
we accept this identification then, for moderate doping, 
Tpair ~ again depending somewhat on the partic- 
ular material being studied. 
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D. Relation to experiments 

1. Estimates of the model parameters 

To begin with, it is necessary to estimate the values of 
the important interactions which determine the behavior 
of the model. The physics is driven by the local anti- 
ferromagnetic correlations between spins, so a priori we 
expect the interactions, other than those within a single 
stripe, to be some fraction of Jaf , which in the high tem- 
perature superconductors is in the range lOOOfC — 2000iir 
pO| . For similar reasons, the bandwidth in the environ- 
ment, W, is expected to be a few times Jaf', numerical 
simulations for the square lattice lead to the estimate 
that the hole bandwidth is approximately 2. 2 Jaf- 
On the other hand, a naive estimate of the bandwidth 
W of the IDEG is given by the bare value, 2t leV, 
although this is certainly reduced substantially due to 
virtual (high energy) single-particle excursions into the 
environment, i.e. leakage of the hole wavefunction into 
the insulating neighborhood of the stripe. 

More detailed estimates may be obtained from exper- 
iment. Since e* /2 is the binding energy of a holon in 
the stripe, we expect that it also determines the tem- 
perature at which stripes begin to lose their integrity, so 
we estimate that e* ^ 2T{^. Thus, e* is certainly re- 
markably small, e* ~ Jaf/2, but still large enough that 
the peculiarities of the small s* limit are avoided. Simi- 
larly, if we identify Tpair with the spin-gap temperature 
deduced from NMR, we can approximately invert the re- 
lation Tpair ~ J^sp/^* to obtain an estimate of Jgp ~ e*, 
where the exact numerical relation between these two 
quantities depends on numerical amplitudes which we 
cannot calculate with any great accuracy. For this range 
of parameters, it also follows that Ai(0) ~ Tpair, con- 
sistent with estimates of the superconducting gap from 
photoemission experiments. Finally, from the magnitude 
of the pseudogap observed in c-axis optical response, we 
estimate that « s*. This implies that the cuprates 
lie in the crossover region between large and small e* 
(regimes B and C described in Sec. VI), which is also 
the region of maximum Tpair, as shown in Fig. ^. We 
feel that these values of e* , Jsp, and Ag are physically 
reasonable. 



2. Does local pairing on stripes provide a consistent 
explanation of the pseudogap behavior of underdoped 
cuprates? 

In the above discussion, we interpreted the experi- 
mentally measured pseudogap behavior in underdoped 
cuprates as superconducting pairing in a large range of 
temperatures above Tc. This behavior was predicted by 
us [|| on the basis of a phenomenological analysis of the 
relation between the superconducting Tc and the mea- 
sured zero temperature superfluid phase stiffness (i.e. the 
zero temperature London penetration depth). It provides 
a very natural explanation of the "spin gap" behavior 



that has been widely observed in planar copper NMR 
measurements in underdoped cuprates |^^. Here, there 
is a peak in l/TiT at a characteristic pairing tempera- 
ture above Tc, below which there is a rapid falloff that 
is quite similar to that observed below Tc in more heav- 
ily doped cuprates. The interpretation of the spin gap 
as a superconducting gap has recently received consid- 
erable support from ARPES experiments which 
find that the magnitude and wave vector dependence of 
the pseudogap above Tc is similar to that of the gap 
seen well below Tc in both underdoped and optimally 
doped materials. The temperature above which this gap 
structure becomes unobservable correlates well with the 
pairing scale deduced from spin gap measurements. Mea- 
surements of the in-plane optical response are also highly 
suggestive of super cond ucting pairing above Tc in under- 
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doped cuprates ||9£ 

This interpretation has been questioned because a 
large fluctuation diamagnetism and conductivity have 



not been observed between Tc and Tpair [101 1- How- 



ever, we believe that the absence of dramatic magnetic 
field effects is readily understood. Well above Tc, the 
superconducting fluctuations are essentially one dimen- 
sional, with little effect of the Josephson coupling be- 
tween stripes. Consequently, an applied magnetic field 
does not drive any significant orbital motion until coher- 
ence develops in two (and ultimately three) dimensional 
patches, close to Tc. We are currently engaged in more 
detailed calculations of these effects, to make this argu- 
ment more quantitative. 

Recently it has been determined that in un- 

derdoped and optimally doped La2-a;Sra;Cu04, there is 
a unique relation between the mean separation between 
stripes {i.e. the half-period of the dynamical incommen- 
surate spin fluctuations) and the superconducting Tc. We 
have previously predicted such a relation |86) as a natural 
consequence of the existence of superconducting fluctua- 
tions on a single stripe and the idea that Tc is determined 
by the Josephson coupling between stripes. 



3. Commensurability and Near Commensurability Effects 

The charge density on the stripes (and hence, the value 
oi kp) is largely determined by the competition between 
the local tendency to phase separation and the long-range 
Coulomb interaction; however, there are commensura- 
bility effects both within the IDEG (which tend to pin 
2kpa = 21: jm where m is the order of the commensura- 
bility) and transverse to the stripes, which tend to pin 
the spacing between stripes at an integer times the lat- 
tice constant |16|. In La2-a;Sr2;Cu04, neutron scattering 
evidence supports the notion that there is a strong ten- 
dency toward locking the hole density within a stripe 
near commensurability m = 4 for a range of x less than 
X = 0.125, and to pin the spacing between stripes near 4 
lattice constants for x > 0.125. (See Sec. VIIB.) Within 
the theory of the IDEG, commensurability leads to a 
charge gap and insulating behavior. However, for a weak 
commensurability, the gap develops at low temperatures 
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where it must compete wi th s uperconductivity. (For an 
alternative view, see Ref. 



log] 



4- Are there any experimentally testable predictions that can 
be made on the basis of this mechanism? 

To begin with, it is important to stress that there al- 
ready exists considerable experimental evidence that the 
physics discussed in this paper is pertinent to the high 
temperature superconductors. Some of this has been dis- 
cussed above. Neutron scattering and transport mea- 
surements provide direct evidence of hole-rich metallic 
stripes in an antifcrromagnetic environment in at least 
the La2Cu04 family of materials. The convincing exper- 
imental evidence that underdoped cuprates behave like 
granular materials in that a superconducting gap opens 
well above Tc, strongly suggests that the superconduc- 
tivity is inhomogeneous at some intermediate scale of 
length and time. Moreover, the absence of strong effects 
of magnetic fields in a regime of strong superconduct- 
ing fluctuations indicates that these inhomogeneities are 
likely to be one-dimensional in character. The fact that 
both s-wave and d-wave symmetry are manifest in differ- 
ent phase-sensitive experiments on essentially the same 
materials supports the idea that there are strong, local 
fluctuations which break the (approximate) four-fold ro- 
tational symmetry of the crystal |103|j87| ] . 

However, while we feel that these experimental facts 
provide strong evidence for the general form of our model, 
they do not probe microscopic structure of the proposed 
pairing mechanism. There are, however, various signa- 
tures that could, in principle, be detected. We predict a 
spin-1, charge zero excitation (a quasi one-dimensional, 
magnon mode) with an energy gap A^, which is of or- 
der the superconducting gap. This mode could, in prin- 
ciple, be detected in neutron scattering. We also pre- 
dict one or more gapped charge=0, spin=0 modes, the 
breathers; for the expected case of a — 2/3 there are 
two such modes, and the lowest energy one should also 
have energy gap As. This mo de co uld, in principle, be 
observed by Raman scattering |104]. Since it also could 
hybridize with a phonon, it could also show up in neu- 
tron scattering. It is interesting to note that a phonon 
mode wh ich is sensitive to the onset of superconductivity, 
[104,105| and a magnon, both with energy about 40meV 
have been observed [106| in the superconducting state of 
YBa2Cu307_5; we are currently exploring whether these 
two phenomena reflect the two collective modes discussed 
above. 

A stripe structure may have a nematic phase, in which 
the stripes are orientationally ordered along a partic- 
ular direction. Such a phase should display a strik- 
ing anisotropy in its phase stiffness. It is interesting 
to note that a big increase in the phase stiffness is ob- 
served as YBa2Cu307_5 is overdoped [107|. This behav- 
ior has been attributed to superconductivity (induced by 
the proximity effect) in the CuO chains, as they become 
filled. However, such an interpretation requires that the 
superfluid density in the chains is greater than in the 



planes, where it originated. Experimentally it may not 
be easy to distinguish nematic stripe order in overdoped 
YBa2Cu307_5 given the existence of the CuO chains. 

One feature of our model is that there are two, physi- 
cally distinct, spin gaps, one associated with the IDEG, 
and hence with the "superconducting gap" , and the other 
(larger gap) with the insulating environment. However, 
in practice, we expect that the two gaps will be similar 
in magnitude because the difference will be "smoothed 
out" by the motion of the holes between the stripe and 
the environment. (Exactly this sort of "smoothing out" 
of the gap occurs in the "Cooper limit" for the conven- 
tional proximity effect.) Finally, we observe that there 
are calculable consequences of our model for single par- 
ticle properties, such as the density of states, which are 
currently under investigation. 

Another qualitative test of our ideas is to look for 
high temperature superconductivity in new materials 
that have one-dimensional metallic and spin-gapped re- 
gions in close electrical contact built into their structure, 
and not necessarily self-organized. In this regard, we note 
that a material with even-leg undoped ladders (which 
have a spin gap ||2^ ) in intimate contact with doped 
Cu02 chains should display the mechanism of supercon- 
ductivity that we have proposed here. Interestingly, su- 
perconductivity with Tc = 12K has been observed |10S] 
at a pressure of 3GPa in Sro.4Ca13.6Cu24O41.84, a mate- 
rial with this kind of structure, although the chains and 
ladders are in different planes, so the electrical conatct is 
not as strong as we would like. At atmospheric pressure, 
it appears that the doped holes are in the chains |10£] 
but, at present, it is not known if this feature persists at 
the high pressures required for superconductivity. 

Our model also could be studied by numerical tech- 
niques. In particular, an environment with a spin gap 
could be represented by either a two-leg ladder or an 
incommensurate dimerized half-filled chain. An envi- 
ronment without a spin gap would be a half-filled one- 
dimensional Hubbard model. In either case the coupling 
to the IDEG should involve strong single-particle or pair 
hopping and a repulsive interaction between holes. 
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APPENDIX A: PERTURBATIVE 
RENORMALIZATION GROUP ANALYSIS 

There are three related senses in which we use the 
renormaUzation group to analyze a complex physical 
problem, such as the present one. 

1) Firstly, the rcnornmalization group, and in partic- 
ular the notion of fixed points, is a theory of theories, 
and it provides a context and structure which allows the 
problem to be approached in the context of its global 
phase diagram. Even when calculations are not carried 
out by use of the renormalization group, the results are 
fundamentally informed by its structure. For instance, 
so long as an exactly solvable model, and a particular 
problem of physical interest are governed by the same 
fixed point, the solvable model can be said to be an ac- 
curate representation of the low-energy physics of the 
problem of physical interest, whether or not there is a 
microscopic correspondence. It is in this sense that a 
large class of physically diverse one-dimensional systems 
can all be described as "Luttinger liquids" , or that the 
resonant level model represents a solution of the antifer- 
romagnetic Kondo problem. Similarly the exact solution 
of the pseudospin model, presented in Sec. V describes 
the physics of the paired spin liquid phase of the IDEG 
in an active environment. 

2) The notion of an unstable fixed point (or line of 
fixed points) also underlies the use of field theories to de- 
scribe condensed matter systems. Of course, condensed 
matter systems have a finite lattice spacing. However, 
in the proximity of an unstable fixed point, the correla- 
tion length diverges, so that the continuum limit is ac- 
tually realized when the correlation length diverges, but 
this is equivalent to holding the correlation length fixed, 
and letting the bandwidth diverge, as is done in defin- 
ing a field theory. Thus, all the field theory results we 
employ, incuding the results based on the equivalence be- 
tween different field theories which goes under the title 
of bosonization, are based on the proximity of the system 
to the Luttinger liquid line of unstable fixed points. 

3) The renormalization group is also a computational 
scheme, which in most cases must be carried out in the 
context of a perturbative evaluation of the beta function. 
The terms "relevant" or "irrelevant" in the renormaliza- 
tion group sense refer to the results of a perturbative 
evaluation of the beta function in the neighborhood of 
a particular fixed point. Such methods are useful for 
determining the stability or lack thereof of a particular 
fixed point. However, in the case in which there is one or 
more relevant interaction, these results can only be used 
to guess the nature of the actual ground state. 



1. Perturbative treatment of Hint 

One approach to the problem is to treat Hint as a 
small perturbation. Thus, one imagines determining the 

properties of the fixed point corresponding to the decou- 
pled problems of the IDEG and the environment, and 



then assessing the relevance of Hint at that fixed point. 
Because, by assumption, the environment has a charge 
gap, any interaction involving excitations of the charge 
degrees of freedom of the environment is irrelevant in 
the renormalization group sense. Thus, Hpair and the 
charge and charge-current interactions in Hint (i.e. the 
terms proportional to Vc and Jc) are immediately seen 
to be irrelevant. In the case in which the environment 
has a preexisting spin-gap, the same analysis implies that 
the remaining interactions in Hint are also perturbatively 
irrelevant. Even in the case in which the environment 
has gapless excitations (^i > 0), the spin couplings can 
readily seen to be perturbatively irrelevant. Thus, for 
weak enough coupling between the IDEG and the envi- 
ronment, the coupling can be ignored in the sense that 
the low energy behavior is qualitatively similar to that of 
the two subsystems in the absence of their coupling. 

In the problem of physical interest, the energy to trans- 
fer a pair of holes from the IDEG to the environment, 
£* , is very small compared to the bandwidth. As we 
have shown in the main body of the paper, this implies 
that the perturbative analysis about the Hint = fixed 
point is valid only in an extremely restricted regime of 
parameter space. In particular, for fixed small, but non- 
vanishing tsp, there is a critical value of e*, such that 
'Hpair is irrelevant for e* > Sc, and relevant for s* < Sc- 



2. Perturbative RG about the non-interacting fixed 
point 

The standard ("g-ology") treatment of the IDEG, may 

be derived by computing the beta function in powers of 
the interactions, ga, using a version of Anderson's poor- 
man's scaling, in which states at the band edge are inte- 
grated out, and new effective interactions are computed 
for the model with a reduced bandwidth, E < W . The 
variation of the coupling constants as a function of E 
are determined by a difi'erential equation, in which the 
microscopoic values of the interactions serve as initial 
conditions. This method can only be applied if all the in- 
teractions are weak on the scale of the bandwidth, as it is 
based on perturbation theory about the non-interacting 
fixed point. 

For the present problem, one can similarly derive the 
appropriate scaling equations for the entire set of inter- 
actions in perturbation theory about the non-interacting 
fixed point. To do this, we notice that the model de- 
fined in Section II is a particular form of an asymmetric 
two-band model, with appropriate couplings, and with 
bandwidths W and W, respectively. However, because 
of the large difference in the bandwidths, the integrating 
out of high energy degrees of freedom, which is the busi- 
ness end of this sort of calculation, must be carried out 
in two stages. In the initial stages of renormalization, 
we integrate out degrees of freedom (of the IDEG) with 
energies between W and E, where W > E ^ W. The 
resulting scaling equations apply so long as all the inter- 
actions remain small [i.e. so long as perturbation theory 
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is adequate) until E reaches the scale of W. For fur- 
ther reduction of the bandwidth, excited states of both 
the environment and the IDEG are being simultaneously 
eliminated. In this way, starting with a set of bare cou- 
pling constants, one obtains a set of renormalizcd cou- 
pling constants at the end of the first stage of renormal- 
ization which serve as initial conditions for the second 
stage flow equations. 



a. The RG flows for W > E 

To begin with, we ignore the differences in band- 
width, so that the model is equivalent to the two -band 
model considered by Varma and Zawadowskii |11C|. This 
allows us to adopt their results (obtained using the usual 
methods) ; translated into the notation of the present pa- 
per, the scaling equations can be written as 
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where v = {vp + vp)/2 is the average Fermi velocity. 



v/vp, P = v/vp, 



Us = Vs 
Uc=Vc 
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Jc 
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and there are three additional scaling equations for ga 
which can be obtained from the equations for ga by plac- 
ing tildes on the gaS and interchanging a and (3. Here, 
we have augmented the original equations of Varma and 
Zawadowskii to include the effects of u mkla pp scattering, 
which was done by Balents and Fisher | 111 | . (We correct 
a factor of two error they made in the scaling equations 
for and g^.) Note that we have adopted the opposite 
sign convention for the beta function than Varma and 
Zawadowskii; here, the dot signifies the derivative with 
respect to £ = log[W/E], which is the negative of their 
variable, log [5*]. 

There are several aspects of these equations that are 
worth noting. In the first place, the scaling equation for 
tsp is the weak-coupling version of the more general Lut- 
tinger liquid result given in Eq. (p^); tgp is perturbatively 



relevant only if [a{3gi—gc)+P{3gi—gc)—Wc] is negative. 
We expect that gc is negative (but possibly small), gc is 
negative and grows in magnitude with renormalization, 
and gi is positive, but typically decreases with renormal- 
ization. Thus, we see that the two ways in which tsp can 
become relevant are through the generation of a large Uc, 
or via spin-gap physics of the environment, in which case 
gi is negative and grows with renormalization. That the 
latter possibility is the more robust is further emphasized 
by the expected large value of /?, which means that the 
term involving gi makes the largest contribution to the 
beta function. In either case, by examining the depen- 
dence of the beta functions of the various other interac- 
tions on tsp, it is clear that once tgp becomes sufficiently 
large, the there is a bootstrap effect which accelerates 
the flows to strong coupling, in that a large tsp makes a 
positive contribution to the beta functions for gc, gc, and 
Uc, and a negative contribution to gi and gi. 



b. The RG flows for W > E :^ W 

We now return to the problem of determining the beta 
function for the initial stages of the elimination of high 
energy degrees of freedom. The scaling equations for the 
regime W > E ^ W can be obtained from the above 
equations by taking the limit vp oo; this has the ef- 
fect of projecting out any intermediate states involving 
the propagator in the environment. The result is the 
scaling equations which govern the initial renormaliza- 
tion process: 
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Most importantly from these equations it is clear that, 
in the initial stages of renormalization, tgp is reduced from 
its microscopic value, although if the interactions in the 
IDEG are not too strong, this reduction may not be too 
severe. There is also an additive negative contribution to 
gi and a positive additive contribution to gc generated 
in this initial stage or renormalization. This is a form 
of asymmetric screening which tends to increase the rel- 
evance of tsp in the final stages of renormalization. How- 
ever, it seems to us unlikely that this latter effect is strong 
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enough to make tgp robustly relevant at low energies in 
the absence of an environmental spin gap. 



APPENDIX B: SYMMETRIES OF THE MODEL 
AND THE COMPOSITE ORDER PARAMETER 

1. Symmetries of the Model 

To begin with, we tabulate the symmetries of the 
Hamiltonian of the IDEG in an active environment, Eqs. 
(I). 

• Parity is a Z(2) symmetry of the system, which 
results in the transformation 

i^2,a{x) Ipi^cri-X), (Bl) 

and the analagous transformation for the environ- 
metal operators. In terms of bosonic variables, 

eaix) 9a{~x) 

Mx) ^ -M-x). (B2) 

where a denotes s or c. Under the action of the 
parity transformation, , pc, and ps are even, and 
P^, jc, and js are odd. 

• Time reversal is a second Z(2) symmetry of the 
system, which results in the transformation 

i^i^ix) 'i^2,i{x), 

i^2,^ix) i^i,i{x), 

i^i,lix) -iip2,T{x), 

tp2,i{x) ^ -itpi,T{x), (B3) 

and the analagous transformation for the environ- 
metal operators. In terms of bosonic variables, 

9s{x) ^ Osix) - VV2, 

Mx) ^ -Mx)- (B4) 

and, of course, i —i. Under the action of the 
time reversal transformation pc and js are even, 
P^, jc, and Ps are odd, P^ transforms as P^ — > 
— exp(i7rTO)Pl^. and the corresponding environ- 
mental operators transform in the same fashion. 

• Spin rotational symmetry is respected entirely 
by the model as originally written, so there is a cor- 
responding SU(2) symmetry of the system, which 
transforms the operators according to 

ipx.a ^^[exp{ij ■ a)]^^^,ijx.a'- (B5) 



and the analagous transformation for the envi- 
ronmetal operators. Manifestly, this transforma- 
tion leaves all the charge, charge current, and sin- 
glet pairing operators invariant, and rotates all 
spin-vectors in the appropriate fashion. Abelian 
bosonization of the model obscures this symmetry, 
which is manifest as a non-trivial relation between 
Kg and 171. Generalizing the original model by 
defining distinct couplings gi^± and y would give 
arbitrary values of Ks and gi (which now should 
be identified with 51, in this case, only the U(l) 
symmetry associated with rotations about the z 
axis remain of the original spin-rotational symme- 
try. The full SU(2) transformation is complicated 
in terms of the bosonic variables, but rotations 
about the z axis correspond to an additive phase 
shift to 9s. 

• Gauge invariance or charge conservation, is man- 
ifest as a global U(l) symmetry of the model (since 
we have not explicitly included the gauge fields) 
which transforms the operators as 

ipx^a exp(i7)V'A,cr, (B6) 

and the analagous transformation for the environ- 
metal operators. In terms of bosonic variables. 




and 0a and (pa are invariant. This transformation 
leaves all the particle conserving operators invari- 
ant, and multiplies all pairing operators by a factor 
of exp[— 2z7]. 

• Translational (chiral) symmetries: There are 
the two independent symmetries corresponding to 
translations (chiral transformations) of the IDEG: 

ipi^a exp(z7t)-0i,o-, 

^2,cr ^ exp{-ijt)tp2,<j, (B8) 

and the analagous tranformations, defined in terms 
of a second, independent angle 74, for the environ- 
mental operators. In the absence of umklapp scat- 
tering, {i.e. if we set gz equal zero) then 74 can take 
on any real value between and 27r, i.e. there is an 
additional U(l) symmetry associated with transla- 
tions of the IDEG). In terms of bosonic variables. 




and the analagous relations (with 7f) for the envi- 
ronmental operators. 

• Spin chiral transformations 

There is an analagous transformation, which 
amounts to a translation of the SDW fluctuations 
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by a half a period, in which the up and down spin 
components are translated in opposite directions. 
We define the spin chiral transformation, C as 

V'2a^#2a, (BIO) 

which in terms of the bosonic variables is, 

^s^cps + (Bll) 

and we define the analagous transformation for the 
environmental operators as C. Hujeg is invariant 
under C, but it has the effect of rotating ps and 
js by TT about the z axis and changing the sign of 
both P''' and Pg, so it is not a symmetry of the 
full Hamiltonian; however CC manifestly is. Hav- 
ing said this, it is clear that additional symmetries 
can be constructed by combining C and C with 
spin rotations by tt about the z axis; we call these 
transformations R and R, and they correspond to 
shifts of 6s and 9s by respectively. In this 

way, an additional discrete group of related sym- 
metry transformations can be constructed consist- 
ing of the identity, CC, CR, CR, CR, and CR; 
this group is Abelian, with a simple group multi- 
pliction table, which is readily obtained. Notice 
that, as with time reversal symmetry, this group's 
operation on spinor fields is double valued. 

• T symmetry There is one additional hidden Z(2) 
symmetry of the Hamiltonian, which combines spin 
and charge transformations, and which is the sym- 
metry that is spontaneously broken in the paired- 
spin-liquid state. This symmetry is combines a 
spin-chiral transformation of the IDEG, C, a tt 
rotation of the environmental spins, P, and an 
inequivalent gauge transformation of the charge 
modes of the IDEG and the environment. In terms 
of the fermionic fields, this symmetry corresponds 
to the transformation 

^2,i ^ *V'2,i- (B12) 

In terms of bosonic variables, this transformation 
takes 





(B13) 



This transformation leaves pc, jc, Pc and jc, in- 
variant, rotates ps, is, P^, and jg by tt about the z 
axis, changes the sign of P^ and P^, and transforms 

'^'^^^ ^m' 

In the above, it is important to realize that a shift 
in the bosonic phases (j)a by ±\/it/2 is equivalent to a 
displacement through a distance equal to the average 
spacing between the particles. For 0c (0s )i spins-cr are 
displaced in the same (opposite) direction. This shift 
leaves the Hamiltonian of the IDEG unchanged because 
the arguments of the cosines in the gi cos('\/87r0s) and 
93 cos(V87r(/)c) terms are changed by 27r. To appreciate 
the significance of this observation, consider the ground- 
state degeneracy of the IDEG with a half-filled band. A 
shift of either ipc or tps by ± ^7r/2 changes the sign of the 
operator o-V'i.o-? since its boson representation is pro- 
portional to exp [i-\/27r(0c + f0s)] . Thus, if this operator 
is ordered the ground state is two-fold degenerate. This 
occurs if both gi and 173 are relevant, as for example in 
the negative-?/ Hubbard model with additional nearest 
neighbor repulsions V , and it is easily understood from 
a strong-coupling analysis, as the ground state is a pe- 
riod 2 charge density wave. These considerations must be 
taken into account in studying the full symmetry group 
of the IDEG as they imply that not all the symmetry 
operations discussed above are linearly independent. 

2. The non-local order parameter 

The non-local order parameter defined in terms of the 
unitary transformation in Eq. ( p9|) is 

Ocomp = UP^U^ (B14) 
= (tto)""^ exp[i\/27r((?c — ^c) cos[\/27r0s] 

can be expressed as a non-local function of the original 
fermionic fields as 

Ocomp = exp [iTT / dyjc{y)]P^- (B15) 

Clearly, this composite order parameter is odd under t 
symmetry. 

APPENDIX C: THE NATURE OF THE "PAIRED 
SPIN LIQUID" 

Various definitions of a "spin liquid" are used in the 
literature ||2^. Here, we define a spin hquid to be a 
quantum disordered ground state of the spin degrees of 
freedom of a system, which means that spin-rotation in- 
variance is unbroken. We also require that translation 
invariance be unbroken for the system to qualify as a 
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liquid. In addition, to distinguish the spin liquid from a 
quantum paramagnet and a Fermi liquid, we require that 
a spin liquid support spinon excitations in its excitation 
spectrum. 

The ground state of a spin-1/2 Heiscnberg chain is a 
gapless spin liquid [ p8| . An integer spin chain or a even- 
leg half-integer spin ladder fail to qualify because spinous 
are confined. (The only finite energy states are integer- 
spin magnons; spinous are bound by a linear potential 
in pairs, or to the ends of chains |112|.) The frustrated 
spin-1/2 chain {e.g. the Majumdar- Ghosh model |24|) 
fails to qualify because translational symmetry is spon- 
taneously broken in the ground state. (See Appendix B.) 
The IDEG away from half-filling displays two kinds of 
behavior; a) when gi is irrelevant, it is a gapless spin liq- 
uid in the universality class of the spin-1/2 Heisenberg 
chain, b) when gi is relevant it has a gap because of 
spinon pairing and is in the universality class of doped 
polyacetylene |113| or a doped Majumdar- Ghosh model 
p]| . It is this latter case, in which spinon pairing causes 
a gap or pseudogap in the spinon spectrum, that we call 
a "paired spin liquid" ; spinous are paired in the same 
way as electrons in a superconductor, and they must 
be created in pairs, i.e. by breaking a bound pair which 
exists in the "vacuum" . 

There are, to the best of our knowledge, only two other 
theoretically well established examples of a spin liquid, 
according to the above definition. The first is the super- 
conducting state of charged particles in high er dimen- 
sions; in this context, it has been shown |114] that the 
usual Bogoliubov quasiparticles have spin 1/2 and charge 
0, where both quantum numbers are sharp quantum ob- 
servables. Clearly, the pairing of spinous in the super- 
conducting state is precisely the pairing that gives rise 
to superconductivity. However, while this connection is 
useful for intuitive purposes, we feel that this state should 
probably not be referred to as a spin liquid, and so we 
propose adding to the above definition of a spin liquid the 
condition that large-scale gauge invariance (in the usual 
sense of superconductivity) should also be an unbroken 
symmetry of the ground state. The second example is 
afforded by some quantum Hall liquid states of electrons 
with spin [115|. For instance, in a quantum Hall system 
consisting of a Laughlin liquid |116| of strongly-paired op- 
posite spin electrons at filling factor v — 2, it is easy to 
see that there exist quasip artic les with spin 1/2, charge 
0, and semionic statistics |117]. This sort of state is a 
realization of the so-called chiral spin liquid 
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